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Preface 


These  notes  are  intended  for  students  with  no  previous  experience 
with  Green’s  functions.  A  knowledge  of  applied  mathematics  at  an  advanced 
undergraduate  level  is  assumed;  for  example,  it  is  assumed  that  the  reader 
will  be  familiar  with  contour  integration  and  with  expansions  of  functions 
of  several  variables  in  the  usual  sets  of  orthogonal  functions.  A  knowledge 
of  integral  transform  methods  would  also  be  useful.  However,  a  brief 
introduction  to  the  Fourier  and  Laplace  transforms  is  given  in  one  of  the 
appendices.  Since  a  wide  variety  of  conventions  is  in  use  this  also  serves 
to  define  the  conventions  and  notations  used  in  the  text. 

This  work  was  sponsored  by  the  Institute  of  Ocean  Science  and 
Engineering  at  The  Catholic  University  of  America  for  use  by  workers  in 
acoustics.  This  interest  is  reflected  in  the  chapter  on  the  Helmholtz 
equation  (Chapter  5)  where  the  applications  are  all  chosen  from  acoustics. 

For  readers  with  an  interest  in  this  field  but  with  no  previous  knowledge 
of  Green's  functions  it  is  suggested  that  the  notes  be  read  from  the 
beginning  with  the  possible  exception  of  the  chapter  on  the  diffusion 
equation  (Chapter  3). 

The  material  in  its  present  form  is  considered  to  be  a  preliminary 
presentation.  It  was  felt  that  it  would  be  desirable  to  make  this  available 
fairly  quickly  and  to  write  a  second  edition  after  some  feedback  had  been 
obtained.  Consequently,  the  author  will  be  especially  grateful  to  any 
readers  who  will  take  the  time  to  offer  criticisms  and  corrections. 

The  author  acknowledges  with  gratitude  the  support  of  the  Office  of 
Naval  Research  (Code  468) . 
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I.  XSTKODEJCTION 

We  will  introduce  Green's  function  by  means  of  a  strap le  example,  and  in 
later  chapters  discuss  some  particular  equations  in  detail.  To  explain  our  choice 
we  first  review  sons  general  properties  of  second  order  linear  partial  differential 
equations. 

1.1  Classification  of  partial  differential  equations 

Any  discussion  of  second  order  equations  oust  begin  with  their 
classification  as  elliptic,  parabolic,  or  hyperbolic,  and  the  prototypes  of  these 
categories,  the  Laplace,  diffusion,  and  wave  equations: 

v2* 
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Mathematically,  the  classification  arises  from  questions  of  existence,  the  problem 
being  to  determine  boundaries  and  boundary  conditions  such  that  existence  theorems 
can  be  established.  We  review  only  the  general  outlines  of  this  scheme  here,  and 
leave  more  precise  statements  for  later  chapters. 

A  condition  which  specifies  the  values  of  a  function  on  a  boundary  is 
called  a  Dirichiet  condition.  If  the  derivative  (  *  V  )  of  the 

function  in  a  direction  n  normal  to  the  boundary  is  given  the  condition  is  called 
a  Neumann  condition.  A  linear  combination  of  the  function  and  its  normal  derivative 
is  called  a  mixed  condition,  while  specifying  these  independently  is  called  a 
Cauchy  condition. 


For  an  elliptic  equation  Dirichiet,  Neumann,  or  mixed  conditions  on  a 


closed  boundary*  are  appropriate.  For  example,  the  steady  state  temperature 
distribution  in  a  solid  satisfies  Laplace’s  equation.  A  Dirichlet  condition  would 
give  the  temperature  on  the  surface,  while  a  Neumann  condition  would  give  the  heat 
flux  through  the  surface.  (By  Fourier’s  law  the  heat  flux  is  given  by  H  =  kVT,  k 

■y 

constant.  Therefore,  the  normal  component  of  H  is  H  «*  kn*7T  »  k  3T/9n.) 

n 

If  a  Neumann  condition  is  given  it  is  evident  that  an  additional  condi¬ 
tion  must  be  imposed,  for  if  the  heat  flux  were  everywhere  into  the  body  it  is 
clear  that  a  steady  state  could  not  be  reached  and  the  boundary  value  problem 
would  have  no  solution.  Therefore,  we  impose  the  additional  condition  that  the  net 
heat  flux  through  the  surface  vanish,  i.e., 


J  d8  1  "  0 

For  a  hyperbolic  equation  an  open  boundary  is  needed.  Cauchy  conditions 
are  usually  appropriate  over  at  least  part  of  the  boundary,  while  Dirichlet, 
Neumann,  or  mixed  conditions  my  be  given  over  the  remainder.  For  example,  the 
displacement  of  a  vibrating  string,  $  say,  satisfies  the  wave  equation. 

If  the  string  is  clamped  at  x  •"  +1 
then  the  boundary  is  an  open  rectangle 
in  the  x-t  plane  with  Dirichlet  con¬ 
ditions  0)  given  on  x  *•  +1,  t  >  0 

and  Cauchy  conditions  (the  initial  values 

)  given  on  t  ®  0,  -1  <.  x  <.  1  —  - 


,  .  W 
°f  *. 


*A  region  with  a  closed  boundary  is  completely  enclosed.  The  boundary  may  be 
closed  at  infinity.  If  the  solution  is  required  to  vanish  at  infinity  this  is 
equivalent  to  imposing  a  Dirichlet  condition  on  a  sphere  and  taking  the  limit  as 
its  radius  becomes  infinite.  Closed  boundary  should  not  be  confused  with  the  term 
closed  region,  by  which  we  mean  that  the  region  includes  its  boundary. 


Parabolic  equations  represent  an  intermediate  case  between  the  elliptic 
and  the  hyperbolic.  A  Dirichlet  condition  over  at  least  part  of  an  open  boundary 


is  usually  appropriate,  with  Dirichlet,  Neumann,  or  mixed  conditions  on  the 
remainder-.  For  example,  the  transient  temperature  distribution  in  a  rod  is 
governed  by  the  diffusion  equation  « 

III 
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0,  a  constant. 

3x“  a* 

If  the  rod  lies  between  x  ■  +  1,  the  boundary  in  the  x«t:  plane  is  an  open  rec¬ 


tangle  with  a  Dirichlet  condition  (initial  temperature)  given  on  t  *  0,  -1  <.  x  <_ 1. 
On  the  remainder  of  the  boundary  x  ■  +  1,  t  >  0  a  Dirichlet  condition  (temper¬ 

ature  at  ends),  a  Neumann  condition  (heat  flux  through  ends),  or  a  combination  of 
these  may  be  given. 

In  addition  to  existence  there  are  two  other  requirements  for  a  well 
posed  problem:  uniqueness  and  stability.  By  stability  we  mean  that  the  solution 
depends  continuously  on  the  boundary  conditions,  i.e.,  a  small  change  in  the 
boundary  conditions  implies  a  small  change  in  the  solution.  The  above  rough  sketch 
of  appropriate  boundary  conditions  has  taken  these  requirements  into  account. 

For  a  discussion  of  the  concept  of  a  well  posed  problem  see  Courant  and  Hilbert 
(1962,  pp.  226-232). 

Equations  may  also  be  characterized  by  their  effect  on  the  boundary  con¬ 
ditions.  The  elliptic  operator  is  a  smoothing  operator: •  any  discontinuities  in 
the  boundary  condition,  as  well  as  its  general  form,  are  smoothed  when  continued 
away  from  the  boundary.  In  contrast  with  this  type  of  behavior,  the  hyperbolic 
operator  has  propagation  properties:  any  discontinuities  and  the  general  form  of 
the  boundary  conditions  persist,  are  propagated,  when  continued  away  from  the 
boundary.  Again,  the  parabolic  equation  is  an  intermediate  case.  In  the  example 

given  above  the  initial  condition  is  smoothed  while  the  boundary  conditions  are 
propagated. 
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Physically,  elliptic  equations  usually  describe  static  phenomena.  For 
example,  electrostatics,  magnetostatics,  steady  state  temperatures  are  governed 
by  Laplace’s  equation.  Hyperbolic  equations  describe  dynamic,  time  dependent, 
phenomena,  e.g.,  electromagnetic  and  acoustic  radiation.  Parabolic  equations  are 
often  very  useful  as  an  approximation,  e.g.,  circuit  theory,  but  their  appearance 
in  physics  is  always  suspect.  In  the  example  given  changes  in  the  boundary  con¬ 
ditions  are  propagated  instantaneously,  which  is  decidedly  unphysical. 

1.2  Introductory  example 

To  introduce  the  Green’s  function  associated  with  a  second  order  partial 
differential  equation  we  begin  with  the  simplest  case,  Poisson's  equation 

V2*  -  -4ttp 


which  is  simply  Laplace’s  equation  with  an  inhomogeneous,  or  source,  terra.  A 
convenient  physical  model  to  have  in  mind  is  the  electrostatic  potential 
$(x,y,z)  which  arises  from  a  source  p(x,y,z)  which  is  a  volume 

distribution  of  charge  (charge  per  unit  volume).  We  imagine  that  a  point  charge 
of  unit  magnitude  (unit  source)  is  at  a  point  ro(xo,yQ,zo)  (source  point) 

and  that  the  potential  is  measured  at  a  point  r(x,y,z)  (field  point). 

From  Coulomb's  law  the  potential  is  just  the  reciprocal  distance  between  the  two 
points  (Gaussian  units  are  being  used).  Written  as  a  function  of  r  and  rQ 
we  call  this  potential  the  Green’s  function 


G(r,ro)  - 


l** 

r-rn 


In  general,  a  Green's  function  is  just  the  response  or  effect  due  to  a  unit  point 

source.  We  also  note  the  symmetry  property  (reciprocity  relation) 

6(r,ro)  -  G(?o,?0 

Suppose  that  there  is  a  charge  distribution  p(r)  in  a  certain  region 

■> 

R  of  space.  The  potential  at  a  point  r  inside  or  outside  of  K  can  be  written 


“4- 


${r) 


P(V 

|r-r  ! 


f  dx  p<r  )  G(r,r  ) 
j  0  0  0 

R 


(dto  =  dx^  dy^  dz^  )  This  results  from  simply  treating  each  element 

of  charge  p(rQ)  as  a  point  charge  and  taking  the  sum  of  the  con- 

tributions  from  all  elements.  The  integrand  is  singular  if  r  is  in  R,  but  the 
integral  ie  convergent.  Thus,  wa  have  obtained  a  solution  of  Poisson’s  equation  by 
means  of  G, 


So  far  this  is  intuitively  obvious  and,  therefore,  unremarkable.  But 
what  comes  next  is  quite  remarkable.  We  now  imagine  that  in  addition  to  the 
source  P  we  also  have  present  certain  boundaries  on  which  boundary  conditions 
are  specified.  Since  Poisson's  equation  is  elliptic,  we  know  that,  for  example, 
a  Dirichlet  condition  is  appropriate.  Therefore,  let  us  suppose  that  the 
equation  holds  in  a  closed  regio.  R  and  that  the  value  of  9  is  specified  as 
a  continuous  but  otherwise  arbitrary  function  on  the  surface  a  of  R.  We  assume 
a  continuous  boundary  condition  for  simplicity  in  this  introductory  example; 
later  we  will  be  interested  in  relaxing  this  condition. 

We  anticipate  two  very  surprising  and  pleasing  results.  The  first  is 
that  the  definition  of  G  (  r,r  ),  given  above  in  the  boundary-free  case,  can  be 
extended  simply  and  used  to  obtain  a  solution  of  the  boundary  value  problem;  the 
second  is  that  the  reciprocity  relation  G  (  r>rQ)  -  G  (  r0>r  )  continues  to  hold 
for  the  Green's  function  so  defined. 

We  have  defined  G  in  the  boundary-free  case  as  the  response  to  a  unit 
point  source.  This  implies  that  the  equation 

V2  G(r,rQ)  -  — 4tt  S(v-rQ)  (1) 

is  satisfied,  and  we  assume  that  it  continues  to  be  satisfied  when  boundaries  are 
present.  This  will  be  part  of  the  definition  of  G.  The  source  term,  our  unit  point 
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source,  is  a  delta  function 


S(r-rQ)  =  6(x-xo)5(y-yo)o(z-zo) 

In  our  electrostatic  model  6(r-v  }  is  the  volume  distribution  PW  of 

o 

a  point  charge  at  .  We  shall  use  delta  functions  freely  as  the  simplest  formal 
device  for  manipulating  Green’s  functions.  An  introduction  to  the  formalism  of 
delta  functions  is  given  in  Appendix  I.  We  recall  the  formal  rules 

S(x-x  )  ■  0  x  r  x 

0  o 

CO 

f(x)  <5(x-x  )dx  ■  f(x  ) 
o  o 

—CO 

where  f  (x)  is  a  continuous  function.  The  latter  is  the  so-called  "sifting  property". 

->  *■>* 

Thus,  we  may  say  loosely  that  6(r-rQ)  is  zero  except  at  r  “  r0  where  it  is 
infinite  in  such  a  way  that 

dr  6(r-r  )  -  1 

o  o 

and  hence,  represents  a  unit  source. 

We  have  emphasized  the  word  formal  because  within  the  context  of  the 
classical  theory  of  real  variables  there  is  no  function  with  these  properties,  A 
function  is  defined  only  if  it  has  a  definite  value  for  each  point  in  its  range 
(infinity  is  not  a  definite  value).  Therefore,  to  make  the  use  of  these 
functions  rigorous  a  broader  mathematical  context  is  needed.  The  best  known  of 
these  is  the  theory  of  distributions.  An  admirable  introduction  to  this  theory 
is  the  book  of  Lighthill  (1958). 

We  shall  not  be  concerned  with  justifying  the  use  of  delta  functions 


because  ail  of  our  results  can  be  obtained  sore  laboriously  without  their  use.* 
Moreover,  ve  regard  the  delta  function  not  only  k  a  convenient  shortcut,  but  also 
a s  supplying  useful  insight  because  of  its  interpretation  as  a  unit  point  source. 
Of  central  importance  in  our  work  will  be  Green's  theorem 


f  dT  -  W2!]}  -  |  d*  (U  f|  -  V  f|} 


which  is  valid  if  U  and  V  are  continuous  with  continuous  partial  derivatives  an  to 


the  second  order  in  the  closed  region  R.  n  is  the  unit  outward  normal  to  R  and 
— 1  ■  n-v  is  the  derivative  in  this  direction.  This  follows  from  the 


divergence  theorem  (Gauss'  theorem) 

f  dr  div  F  | 


n*F  as 


by  substituting 


and  expanding 


div  (U7V) 


VU*W  + 


by  the  well  known  vector  identity,  interchanging  U  and  V,  and  subtracting. 

Returning  now  to  our  boundary  value  problem  we  write 

7  <{>  «  -47?  p(r  ) 

0  o 


where 


7q  G(r,rQ)  »  -4n  5(r-rQ) 


2  2  2 

— -  +  2  +  ~ 

3x  3y  3z 

O  0  o 


*&$  an  example  of  this  approach  we  cite  the  book  of  Soamerfeld  (l9^9)j  concise, 
rwjcrous  proofs  of  the  theorems  needed  are  given  by  Titchmarsh  (1958,  chap.  21). 
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In  the  first  equation  ve  have  simply  replaced  r  by  rQ.  The  second  follows  from 
eq.  (l),  the  symmetry  of  the  Green's  function,  and  the  fact  that  the  delta  function 
is  even, 

6(x-x)  -  S(x-x) 

o  o 

However,  since  ve  have  not  fully  defined,  let  alone  proved  the  symmetry,  of  Green's 
function,  ve  night  do  better  to  simply  interchange  tbs  roles  of  r  and  tq  in  our 
interpretation  of  the  G  in  eq.  (3)  and  regard  it,  as  the  response  at  to  a  source 
at  r  ,  Really,  an  interpretation  at  this  point  is  not  necessary,  since  what 
follows  is  sinpiy  fores!  manipulation.  We  could  just  as  easily  take  eq.  (3)  as 
part  of  the  definition  of  G  and  later  prove  the  reciprocity  relation  and 
derive  eq,  (l). 

Multiplying  eq.  (?)  by  G,  eq.  (3}  by  $  ,  integrating  and  subtracting 

vs  obtain 

[  dT  {GV2£  -  $V2G  +  An  Gp(?  )  -  At -<f  6(r-r  )}  »  0 

J  o  o  o  o  0 

R 

By  Green's  theorem 

f  dv  {GV2$  -  $V2G)  -  I  ds  (G  ~  -  «> 

J  o  o  o  I  o  Sn  Sn 

J  *  Q  O 

R  a 

and  by  the  sifting  property 

dt  $(r  )  S(r-r  )  «  $(r) 

00  o 

Therefore,  we  have 

»(J)  -  ]  p(r  )«(?.?  )  +  ^  j  dS(j  {G  |f-  -  4  |f-  }  (« 

J  J  o  o 
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wh-’ch  expresses  $  in  terms  of  its  values  and  the  values  of  its  normal 

derivative  on  the  surface.  This  is  not  a  solution  of  the  boundary  value  problem 

because  both  and  3$/ 3n  are  not  given  on  the  surface.  However,  if  we  can 

impose  a  further  condition  on  G,  viz.,  that  it  vanish  on  cf,  we  can  eliminate  the 

term  involving  3<J>/3n  .  This  additional  condition  completes  the  definition 

o 

of  G,  We  would  then  have 

*(r)  -  f  di  p(?o)G(?,?o)  -  ^  |  daQ  <K?Q)  ff-  (5) 

i  •  o 

R  a 

which  expresses  <{>  in  terms  of  its  boundary  values.  We  have  obtained  this 
relation  assuming  that  $  is  a  solution  of  Poisson’s  equation.  We  may  now  ask 
if  we  substitute  for  $  in  the  surface  integral  an  arbitrary  continuous  function 
does  this  relation  give  a  solution  of  the  boundary  value  problem?  It  does  if  a 
solution  exists,  for  if  it  exists  it  must  satisfy  this  relation.  It  is 
Important  for  us  to  be  aware  of  the  appropriate  existence  theorems.  We  may 
conclude  that  ea.  (5 )  gives  a  solution  of  the  Dirichlet  problem  because  there 
is  an  existence  theorem  (cf.chap.  2),  but  eq.  (4)  does  not  give  a  solution  of 
the  Cauchy  problem  for  Poisson’s  equation  because,  in  general,  a  solution  does 
not  exist. 

We  have  reduced  the  original  problem  to  the  simpler  problem  of  finding 
the  Green’s  function,  i.e.,  of  finding  the  potential  of  a  point  change  subject 
to  the  condition  of  zero  potential  on  the  boundaries  (grounded  boundaries). 

There  are  two  general  methods  of  constructing  Green’s  functions,  which  we  will 
explore  in  the  next  chapter.  They  are  the  method  of  expansion  in  eigenfunctions 
and  the  method  of  reflection  or  imaging.  The  latter  method  works  only  for  some 
rather  special  geometries,  but  when  it  is  available  it  provides  a  simple 
analytical  expression  for  the  solution,  with  obvious  advantages  over  an  infinite 


series. 


Having  interpreted  G  as  the  potential  of  a  point  charge  in  the  presence 

of  grounded  conductors,  we  can  illustrate  the  reciprocity  relation  G(r,rQ)  -  G(ro,r) 

as  follows.  Imagine  that  a  point  change  is  placed  at  rQ  near  a  grounded 

conductor  and  that  the  potential  is  measured  at  r  .  This  potential  is  due  to  the 

point  charge  and  a  surface  distribution  of  charge  on  the  conductor  which  is 

brought  up  from  the  earth  and  distributed  so  as  to  maintain  the  potential  of 

the  conductor  at  zero.  If  now  the  point  charge  is  moved  to  r  and  the 

■> 

potential  measured  at  rQ  ,  the  surface  charge  on  the  conductor  will  be 
rearranged  so  that  the  two  measurements  are  identical. 

To  summarize,  the  solution  of  Dirichlet's  problem  for  Poisson's 
equation  is  given  by 

Hi)  -  { dTo  ni0)o(i,i0)  -  ^ }  dso  nij  f|- 

i  o 

R  cf 

where  the  Green's  function  is  defined  by 

VZG  -  -4it  S(r-r  ) 
o 

with  the  condition  that  G  vanish  on  the  boundary.  We  emphasize  that  a  Green’s 
function  is  defined  not  only  with  respect  to  an  equation  and  its  boundary 
conditions,  but  also  with  respect  to  a  particular  region. 
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II.  LAPLACE'S  EQUATION 


The  theory  of  Laplace's  equation  is  usually  called  potential  theory.  The 
classic  work  is  Kellogg's  "Foundations  of  Potential  Theory",  where  rigorous  proofs  of 
the  following  existence  theorems  may  be  found  (chap.  11 ). 

2.1  Existence  theorems 

By  an  interior  problem  we  shall  mean  the  problem  of  finding  a  solution  of 
2 

Laplace's  equation,  7 «  0  ,  within  a  closed  finite  region  R  (the  region 

includes  its  surface  a  )  which  satisfies  boundary  conditions  on  o.  By  an  exterior 
problem  we  mean  the  problem  of  finding  a  solution  in  the  infinite  region  outside 
of  a  surface  subject  to  the  additional  condition  that  vanish  at  least  as  fast 
as  r-1  as  r  -*•  ®  .  This  condition  is  needed  for  uniqueness.  Physically,  it  means 
that  at  sufficiently  large  distances  the  source  distribution  should  look  liice  a 
point  source.  We  have  the  following  theorems. 

There  are  unique  solutions  of  the  interior  and  exterior  Dirichlet 
problems  and  of  the  exterior  Neumann  problem  for  continuous i boundary* values. 

There  is  a  solution  of  the  interior  Neumann  problem,  unique  to  within 
an  additive  constant,  for  continuous  boundary  values  provided 

ds-~£  «  0  . 
a 

To  give  rigorous  proofs  of  these  theorems  one  must  impose  certain 
conditions  on  the  shape  of  the  region  R.  For  a  thorough  discussion  Kellogg's 
book  may  be  consulted.  The  conditions  are  very  weak  from  the  point  of  view  of 
physical  applications  and  we  may  assume  that  physical  boundaries  are  always 
sufficiently  regular. 


11  - 


The  theorens  apply  to  Poisson’s  equation  if  the  source  tern  is  piecevise 
differentiable.  (This  is  not  the  weakest  possible  condition,  but  it  is  a 
convenient  one.) 

2.2  Method  of  reflection 

We  consider  the  Dirichlet  problem  first.  We  have  seen  that  there  is 
a  solution 

$  ■  [  dt  p(r  )G(r,r  )  -  -r^  j*  ds  $(r  )  ■— i— 

J  00*0  4*  j  o  v  o  3n 

R  a 

of  Poisson’s  equation  provided  we  con  find  a  Green’s  function  which  satisfies 

"  -4u  6(r-r^) 


and  vanishes  on  the  boundary.  For  exterior  problems  we  add  the  condition  that  G 
-1 

vanish  as  -  for  large  r  in  order  to  satisfy  the  conditions  of  the  existence 

.-#•  -*■  |  — 1 

theorem.  A  particular  solution  is  lr_r0l  to  which  we  may  add  any  solution  of 
the  homogeneous  equation,  F  say, 

„  1 


Is-  "* 
r-r 


-i-  F(r) 


and  hope  to  choose  F  so  that  the  boundary  condition  is  satisfied. 

The  method  of  reflection  consists  of  choosing  for  F  solutions  of  Poisson 

equation  corresponding  to  point  sources  at  positions  outside  of  the  region  of 

2 

interest  R.  Within  R  these  are  solutions  of  the  homogeneous  equation,  7  G  «=  0 
For  example,  suppose  that  is  given  on  the  plane  z  =  c  and  we  wish 
a  solution  for  the  half- space  z  >0. 

G  can  be  written  as  the  potential  of 

.y 

a  positive  point  change  at  iQ  ( x0>yo>z0 ) 
plus  the  potential  of  a  negative  point 


„ - 

‘  -b 

> 

change  at  an  "image  point"  (xo,yo>-zo)  obtained  by  reflecting  rQ 


in  tee  plane 


G<r‘r0) 


Ir-r 


Ir-rJ 


2.-1/2 


«*  { (x-xq) A  +  (*"Zo)  } 

-  {(x-xq)2  +  (y-yQ)2  +  (2**0)2)  ^2 

The  boundary  condition  is  obviously  satisfied  on  z  =  o. 

This  trick  is  veil  knovn  in  electrostatics  end  can  be  used  for  other 
regions  bounded  by  planes.  For  example,  consider  the  quarter-space  x  >  0,  y  >  0 
We  place  image  sources  outside  this  region  at 
(->  ?1(-xo,yotzo). 


(+)  r2^~x0’-yo,Zo^’  and 

(-)  VvvV- 

The  signs  of  the  sources  are  indicated.  It  is  easily  verified  that  the  boundary 
condition  is  satisfied  on  both  of  the  boundary  planes. 

We  can  think  of  this  as  a  process  of  successive  imaging,  each 
reflection  involving  a  sign  chenge,  until  ve  return  to  rQ  .  As  a  further 
example  (Somznerfeld,  p.8o)  consider  a  60°  wedge. 

There  are  five  imeges,  formed  by 
successive  reflections  as  shown 

in  the  figure.  Obviously,  the  ~*(4-) 

method  is  useful  only  if  there  is  a  finite  number  of  images,  i.e.,  if  ve  do 
eventually  return  to  the  starting  point  rQ  .  For  a  wedge  this  will  be  so  if 
the  wedge  angle  is  tt  times  a  rational  number. 

2.3  Neumanns  problem 

An  obvious  starting  point  is  eq.  (1.4) 

♦©  •  f  +  J  tG  It  -  ♦  !r} 

'00 


Our  first  thought  would  be  to  proceed  as  we  did  for  Dirichlet’s  problem  and 
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Impose  a  condition  on  G,  this  tine  requiring  that  3G/3n  *  0  on  the  boundary. 
However,  for  the  interior  Neumann  problem  this  is  not  possible.  We  can  see  this 
as  follows. 

In  Green's  theorem  let  U  =  1,  V  =  G. 
ds  -|~  *  |  dx  V^G  *  -4*  |  dx  6(r-rQ)  “  — 4-n 

Obviously,  3G/3n  cannot  vanish  everywhere  on  the  surface.  However,  we  can  set 
3G/3n  ■  const.  ■  -4it/I 


on  the  surface,  where  £  is  ^he  total  surface  area.  Then  from  eq.  (1.4)  we 
have  the  solution 


$(r)  -  j  dx  p (r)G(r ,r)  + 

jo  o  o 


4?f  ^ 


G  t —  +  const, 
o  an 

o 


which  reflects  the  fact  that  uniqueness  is  proved  only  to  within  an  additive 
constant  for  the  interior  Neumann  problem. 

The  method  of  reflection  in  a  plane  works  for  Neumann's  problem  if 
it  works  for  Dirichlet’s  problem.  The  only  difference  is  that  there  is  no 
sign  change  under  reflection. 

2.4  Reflection  in  a  sphere 

For  a  spherical  boundary  there  is  an  imaging  principle  for  Diriehlet's 
problem  but  not  for  Neumann's.  The 

image  is  on  the  same  radius  as 

+  2 
the  unit  source  *0  .  Let  ■  a 

where  a  is  the  radius  of  the  sphere  and  let  the  strength  of  the  image  source 


be  -a/r  . 

D 

Then 


a/r 


«-*  -* 
r-r 


|r  -  #->2?  | 
i  'r  o' 
o 
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which  is  valid  inside  or  outside  of  the  sphere. 

2.5  Expansion  in  orthogonal  functions 

The  method  of  reflection  may  produce  an  infinite  number  of  images.  As 
an  example,  suppose  Laplace's  equation  is  satisfied  in  the  rectangle  bounded  by 


x  =  o,a  mid  y-  0,b. 

i 

* 

t  - 

1  ^ 

1 

-  * 

V* 

—  1- 

To  satisfy  the  condition  on 

'  * 

1 

1 

_ l_ 

* 

. *1 

-24. 

Jui. 

x=o  we  reflect  in  this  line. 

-4* 

— 1  _ 

* 

* 

1 _  *1 

i 

A  reflection  in  x=a  satisfies  the 

JL* 

— 1  - 

1  *1 

1 - r 

condition  on  this  line  but  now  the  condition  on  x=o  is  no  longer  satisfied  so  the 
second  image  is  reflected  in  x<=o,  etc.  In  this  way  we  completely  fill  the  plane 
with  images.  For  detailed  examples  of  this  approach  see  Courant  and  Hilbert 
(1953,  PP-  378-386). 

A  more  reasonable  approach  is  to  satisfy  the  boundary  conditions  at 
the  outset  by  attempting  an  expansion  in  orthogonal  functions  each  of  which 
satisfies  the  conditions.  These  are  obtained  by  the  separation  of  variables 
method.  For  the  present  example  they  are  sines  and  the  expansion 

G  -  I  I  A  8ln  22*  sin  J52X 
L.  **,  mn  a  b 

n«l 

is  a  double  Fourier  sine  series. 

Substituting  into  the  equation 

V2G  -  -2tt  6(x-x  )  6  (y-y  ) 

0  0 


which  we  adopt  for  G  in  two  dimensions,  and  using  the  orthogonality  relation 

a 


sin 


mrx 

a 


.  m  vx  , 
sin  — —  dx 
a 


a  x 
2  mo' 


0 
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and  the  sifting  property  to  determine  the  coefficients  km  we  obtain 


Ttab 


m,n 


sin  mux/a  sin  ra7txo/a  ain  tufy/b  sin  ni?yQ/b 
(ffi/a)2  +  (n/b)2 


(1) 


Since  we  know  that  G  is  singular  at  x  =  xq,  y  =  yQ,  we  are  inclined  to 
ask  in  what  sense  it  can  be  represented  by  an  ordinary  Fourier  series.  Justification 
can  be  found  in  the  theory  of  distributions  or  generalized  functions  (Lighthill, 
1958),  We  only  remark  that  as  the  delta  function  has  a  well  known  representation 
as  a  Fourier  integral 


6(x-x  ) 
o 


JL 

2ir 


dk  e 


lk(x-xQ) 


so  a  row  of  delta  functions  (an  infinite  array  of  sources )  can  be  represented  by 
a  Fourier  series,  e.g., 

“  _  HTTX 

l  (-l)m  6 (x-x  -ma)  -  -  £  sin  3 


m*-oo 


n-0 


.  rorx 

-  sin - 

a  a 


(2) 


(cf.  Lighthill,  chap,  5) 

It  may  also  occur  to  the  reader  to  wonder  how  the  expanded  G  when 
substituted  into* 

*  -  -  2?  f  dso  *  If-  (3 

J  0 

can  give  a  solution  which  satisfies  the  boundary  condition  when  each  term  in  the 
expansion  vanishes  on  the  boundary. 


*We  have  adopted  the  convention  of  writing  Poisson's  equation  in  two  dimensions  in 
the  form  -  ~2ttp  .  Therefore,  in  this  formula  we  have  a  factor  of  instead 
of  4tt  and  a  line  integral  in  place  of  the  surface  integral. 


-  16  - 


The  answer  is  that  the  boundary  must  be  approached  from  inside  the  region.  To 
take  a  simple  example,  suppose  the  potential  isf  (x)  on  y  =  o  0  <  x  <  a  and  zero 
on  the  rest  of  the  rectangle.  Then  eq.  (3)  becomes 

a 

*  -  27  j  dXo  f(Xo>  -0  (4) 

0  00 


(the  outward  normal  is  in  the  -  y  direction). 


Prom  eq.  (l), 


ri£_]  »  8  y  sln  SH  sin  y  jL§*SL ™l!k 

3y  Jy“0  a  ^  a  a  ^  2.  ,nbN2 

Jo  Jo  m  n  n  4-  ( — ) 

a 


(5) 


It  can  be  verified  that  the  expansion 


r,  sinh  ktr(l-y/b)  r  n  sin  mry/b 

2  -STH1-  ■  I 

is  valid  for  0  <  y  <  b.  Hence,  as  y  0+  the  inner  sum  in  eq,  (5)  approaches 
ff/2  .  Therefore, 

,  mux 

,3G  ,  4u  v  mux  ,  o 

[t — ]  n  a  —  )  sin - sin  — — 

8y  y  =0  a  -  a  a 

°  m 

y*-0+ 

«*  2u  5 (x-xq) 

from  eq.  (2).  Substituting  into  eq.  (4)  we  see  that  the  boundary  condition  is 
satisfied. 


$  (x,y-KH-) 


dx  f(x  )<$(x-x  )  »  f(x) 

OO  0 


2.6  Discontinuous  boundary  conditions 

Our  formulas  have  been  derived  under  the  assumption  that  $  is  a 
solution  in  a  closed  region  R  (the  region  includes  its  boundary).  This  implies 
continuous  values  of  <S>  on  the  boundary  and  this  condition  also  appears  in  the 
existence  theorems  quoted.  However,  we  notice  that  the  surface  integrals  which 
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occur  in  the  formulas  may  exist  even  when  the  boundary  conditions  are  discontinuous, 
and  so  it  is  natural  to  ask  if  they  can,  in  some  sense,  represent  a  solution. 

It  is  clear  in  this  case  that  there  cannot  be  a  solution  in  the 
closed  region.  However,  we  can  re-define  the  boundary  value  problem  so  that  we 
seek  a  solution  in  the  open  region  which  approaches  the  (possibly  discontinuous) 
boundary  values  as  the  surface  is  approached,  and  an  existence  theorem  can  be 
proved  for  piecewise  continuous  boundary  values  (Tychcnov  and  Samar  ski, 
p.  261  ff.). 

This  may  actually  be  a  more  reasonable  way  to  state  the  problem 
for  certain  physical  applications.  For  example,  across  a  surface  charge 
distribution  u  there  is  a  .Jump  condition  on  the  normal  derivative  of  the 
potential 

.3$.  ,3<K  . 

Wi '  W2  ‘  • 


the  subscripts  indicating  that  the  surface  is  approached  from  one  side  or  the 
other.  3$/9n  is  not  defined  on  the  surface  itself. 

To  justify  our  derivations,  the  conditions  imposed  on  surface  values 
of  functions  for  the  validity  of  Green’s  theorem  (or  Gauss’  theorem,  on  which  it 
is  based)  must  be  weakened.  This  can  be  done;  see  Kellogg  (p.  119 ). 

YTe  now  give  an  example  to  show  how  the  boundary  condition  is  approached 
at  a  point  of  discontinuity.  Suppose  that  Laplace’s  equation  is  satisfied  in  the 

upper  half-plane  y  >  o  and  that  on  y  =  0  we  impose 

'0  -«®  <  x  <  0 

$(x,0)  "  I  0<x<l 

.0  1  <  x  <  « 

We  first  construct  Green’s  function  for  this  region  by  the  method 


of  reflection.  Since  this  is  a  two  dimensional  problem,  the  potential  of  a 
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r.&nge  is  -  log  ]  r-r^  1 


'.ere  f ' 


-log  Ju- 


xo)2+(y-y0)2  + 


log  f(x-> 


x0>^(^7 


[i^-] 

iac  boundary 


SG 


l3y  Jy  -C 

y  A  y  A 


O  0 

(x-x  )“  +  y“ 


{"&  is  the  outward  r.ornal) 
‘  o 


|  dx  $ 
2?  j  o 


1£-  .  If 

3n  r  j 

n  » 


dx 


,  n2  ,  2 

(x-xq)  +  y 


Let  us  consider  the  approach  to  the  discontinuity  at  the  origin  along 


ray  y  =  a  x 
asking  ars  angle  0 
with  the  x  -  axis, 
-  -  tan  6 


The  integral  can  he  performed  giving 


$ 


x  -xll 


Taking  the  limit  x  ■+•  o 


$  -*•  —  {I  +  tan  1  cot9}  ■ 

it  2.  II 


We  see  the  discontinuity  is  not  continued  away  from  the  boundary, 
but  is  immediately  smoothed,  and  that  all  values  in  the  range  of  the  jump  are 
approached  as  0  varies  over  l80°. 

The  approach  of  this  and  the  previous  section  can  be  used  to  justify 
a  method  which  is  given  in  many  textbooks.  Suppose  we  wish  a  solution  of 
Laplace’s  equation  In  a  rectangle  which  satisfies  a  Dirichlet  condition  on  the 
. .undary.  We  break  the  problem  up  into  four  parts  in  each  of  which  a  condition 
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is  given  on  one  side  of  the  rectangle  with  zero  on  the  other  sides.  The 
separate  problems  can  be  solved  by  the  separation  of  variables  method  and  it  is 
claimed  that  the  solution  of  the  original  problem  is  just  the  sum  of  the  separate 
solutions. 

However,  the  separate  solutions  are  given  as  infinite  series  whose 
uniform  convergence  can  only  be  established  if  the  boundary  condition  is 
continuous,  ana,  therefore,  zero  at  the  comers.  But  we  know  from  the  existence 
theorems  that  the  boundary  condition  for  the  original  problem  need  not  vanish 
at  the  corners. 


We  can  justify  a  discontinuous  boundary  condition  at  the  corners 
for  the  separate  problems  by  means  of  the  expanded  Green’s  function  of  52.5. 


If  0(x,0)  -  f  (x)  0  <_  x  <_  a  ,  for  exeagple,  with  $  -  0  on  the  rest  of  the 


boundary,  it  can  be  shown  that  at  the  origin 
<•  -*•  f(0)  along  A 

4  0  along  B 


and  that  $  approaches  intermediate  values  continuously  along  rays  between  A 


and  B. 


The  method  can  he  used  for  other  regions  with  comers,  such  as  boxes 
and  cylinders. 
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2.7  Exercises 


1. 


2. 


3- 


4*  • 


5. 


6. 


7. 


8. 


9- 


10. 


11. 


Prove  the  reciprocity  relation  for  Dirichlet' s  problem.  Hint:  Multiply 
the  equation  for  G(r,r^)  by  GCr,^)  and  vice  versa.  Subtract  the 
equations,  integrate,  and  use  Green’s  theorem. 

Show  that  for  Neumann’s  problem  the  reciprocity  relation  rosy  always  be 
imposed  as  a  separate  condition. 

Find  the  solution  of  Laplace's  equation,  $xx  +  $  **  0,  in  the  half-space 

y  0  which  satisfies  4>^(x,0)  »  g(x).  (Subscripts  indicate  partial 
derivatives. ) 

Replace  the  boundary  condition  of  #3  by  <t>(x,0)  -  f(x)  . 

Find  the  solution  in  the  quarter- space  x  >_  0t  y  >,  0  which  satisfies 


cj>(x,0)  **  0  4>x(0,y)  -  g(y). 

Replace  the  boundary  conditions  of  #5  by 


<Kx,0)  *»  f(x)  $x(0,y)  -  0 

Find  the  solution  of  Laplace's  equation,  <i>  +  <(>  +  * 

“  xx  yy  zz 

space  z  >_  0  which  satisfies  <j»(x,y,0)  »  f(x,y)  . 

Replace  the  boundary  condition  of  #7  by  (x,y,0)  -  g(x,y) 
Find  the  solution  in  the  quarter- space  x  >_  0,  y  il  _a>< 
satisfies  <H0,y,z)  "0  z  >_  0  ^(x,0,z)  »  f(x,z)  x  >_  0. 

Replace  the  boundary  conditions  of  49  by 

<f>  (0,y,z)  =  g(y,z)  y  >_  0  (p (x,0, z)  «  f(x,z)  x  >_  0 

An  infinite  plane  has  a  hemispherical 
boss  of  radius  a.  Find  Green's 
function  for  the  half- space  for  which  this 
surface  is  convex  (see  diagram).  Assume 
that  Dirichlet  conditions  are  given  on 
the  surface. 


0,  in  Che  half- 


z  <  »  which 


Find 
Green's 
Function 
for  this 
region. 
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12.  Assuming  that  Dirichlet  conditions  are  given  find  Green’s  function  for 
the  region  between  concentric  spheres  of  radii  b  >  a.  Hint:  Write 


1  -*• 

+  P(r) 


Expand  F  in  spherical  harmonics  and  use  the  well  known  expansion 

i 

Trr-"  I  wjiBr  TS 

lr-v  r> 

where  r<  (r>)  is  the  smaller  (larger)  of  r,^  and  (0 ,<j>)  *  P^(cos8)e'*'in^ . 


III.  THE  DIFFUSION  EQUATION 

We  must  now  have  a  unit  point  source  in  time,  as  well  as  in  space. 
Therefore,  we  write 

{v2  -hj?  G(^^o;t’to)  "  -*1T  «<r-r0)6(t-t0> 

a 

We  shall  see  that  G  =  0  t  <  tQ  holds  as  well  as  the  reciprocity  relation: 
G(r,ro;t,to)  =  G(?o,r;-to>-t) 

We  emphasize  that  these  are  not  imposed  conditions,  but  are  derived.  They 

are  usually  described  by  the  word  causality,  from  the  interpretation  of  G 

as  the  response  or  effect  due  to  a  cause  represented  by  the  source  term  (a 

pulse  at  t  »  t  ).  Thus,  the  first  relation  says  that  the  effect 
c 

cannot  precede  the  cause,  and  the  second  shows  the  following  symmetry 
between  cause  and  effect  with  respect  to  the  time  scale 


cause 

effect 

cause 

effect 

1 

-t 

-t 

\ 

t 

1 

t 

It  is  quite  reasonable  to  make  this  interpretation  because  time's 
arrow  is  built  into  the  diffusion  equation:  it  is  not  invariant  'under  time 
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reversal.  This  car.  also  be  understood  in  terms  of  stability,  which  we 
mentioned  in  the  introduction  as  one  of  the  requirements  for  a  well  posed 
problem.  The  problem  of  continuing  an  initial  condition  into  the  future  is 
stable,  whereas  continuation  into  the  past  is  unstable.  We  shall  see  that  the 
causal  interpretation  is  somewhat  more  arbitrary  when  applied  to  the  wave 
equation,  which  is  invariant  under  time  reversal. 

3.1  The  Boundary  value  problem 
We  begin  by  writing 


„2,  1  3'1  ,  .■*  „  % 

V  1  at  “  "4lt 

a  o 


V2G  +  «  -4tt  <S(r-r  )6(t-t  ) 

o  2  3 1  o  o 

a  o 


* 


The  second  equation  follows  from  the  reciprocity  relation.  The  operator  on  the 
left  hand  side  is  the  adjoint  of  the  diffusion  operator.  (The  Laplace  and 
wave  equations  are  self  adjoint.)  We  multiply  the  first  equation  by  G,  the 
second  by  <l>  ,  integrate  and  subtract,  giving 


t+ 


dt 


dx  {GV2i|)  -  i|>72G - ^  (il>  ~~  +  G  )  +  4irp  G} 

O  O  O  i.  dt  o  C 

a  o  o 


4it  >p(r,t) 


where  t  +  indicates  the  limit  as  e  -*•  0+-  of  t  +  e  .  We  write 

^  lr +  G  H"  "  tt  (*g) 

0  o  o 

perform  the  t6  integration  for  this  term, 

[G^]^+  =>  - [Giji]  aQ  since  G(t,t+)  *  0 

o 

and  using  Green’s  theorem  obtain 
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<Kr,t) 


F  ato  i 


dr  p(r  ,t  )G 
o  o  o 


l*v 


C 

1  f  dto  I  dso  {G 


v 

0 


3n 


3n 


♦“H  Jd'o“*>t_.0 


4s&  4  o 

We  impose  the  vanishing  of  G  on  the  boundary  if  Birichlet  conditions 

are  given,  and  of  3G/3n  for  Neumam conditions.  Thus,  the  first  term  repre¬ 
sents  the  contribution  from  the  source,  the  second  from  the  boundary  conditions, 
and  the  third  from  the  initial  condition. 

3.2  Boundary- free  case 

We  consider  first  the  case  in  which  there  are  no  boundaries  and 
no  initial  condition.  Then,  having  obtained  this  solution,  we  can  construct 
other  Green’s  functions  using  the  method  of  reflection.  For  Laplace’s 
equation  we  were  able  to  write  this  solution  immediately  be  appealing  to 
Coulomb’s  law,  but  here  the  answer  is  not  so  obvious.  Let  us  begin  with 
the  one  dimensional  equation 


3  G  _  1_  3G 

a  2  "  2  3t 

3x  a 


— 4tt  5(x-x  )6(t-t  ) 
o  o 


Taking  a  Fourier  K- transform  on  x  and  a  Fourier  o>- transform  on  t, 
the  equation  becomes  (see  Appendix  II ) 


<-k2  +  1X, 

a 


ikx  itoi 


-4tt  e 


,  which  is  a  function  of  the  transform  variables  k  and  as,  is  the 


transform  of  G,  i.e., 


w  oo 

jdx  | 


dt  G(x,x  ;t,t  ) 
o  o 
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Solving  the  transformed  equation 


Ikx  tut  ikx  iut 

o  o  „  o  o 

,  n  c  .  2,  e  e 

-uf  v . « —  ”  4ra  i - tt-tt- 

*■^0  iw  ^_2  ,  ,  2t_2 


-  k 


w  +  ia  k 


Inverting  the  ^-transform 


.  “  ikx  -i<o(t-t  ) 

.  4tto  i  [  ,  e  °  e  ° 

V'-V  ■  sr  J  ** - r~T2 — 

'  w  t  ia  k 


The  integranal  has  a  simple  pole  in  the 
lover  half  of  the  complex  “-plane.  We 
may  complete  the  contour  in  the  upper 
(lover)  half- plane  for  t-tQ  >  0 
The  integral  is  zero  for  t-t  <  0  since  no  singularity  Is  enclosed.  For 
t-tQ  >0  ve  have  -2ri  times  the  residue  at  the  pole  (the  path  of  integration 
is  in  the  negative  sense): 


«*  4ira  e 


ikx  -«2k2(t-t  ) 

2  O  0 

e 


(t  >  to> 


Finally,  inverting  the  k-transform 


4rca 

2u 


-ik(x-x  )  -a2k2(t-t  ) 

dk  e  0  e  0 


2a  dk  cos  k(x-x  )e 
o 


-o2k2(t-to) 


2— exp  -  (x-x  )2/4a2(t-0 

SET  0  0 


(t  >  to) 
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G  - 


We  see  that  at  any  instant  after  to  the  effect  of  the  pulse  is  felt  everywhere, 
propagation  is  instantaneous. 

The  method  is  trivially  extended  to  n  dimensions.  We  simply  take  Fourier 
transforms  on  all  variables  and  obtsin 

G  *  4iraZ  (2aA(t-t  ))  n  U(t-t  )  exp  -  |r-r  jV4a  (t-t  ) 

o  o  c  o 


The  unit  function  U,  which  is  equal  to  one  (zero)  for  positive  (negative) 
argument,  appears  because  G  vanishes  for  t  <  t0  .  In  the  older  literature  these 
are  called  the  heat  pole  solutions. 

3-3  Method  of  reflection 

We  are  now  in  a  position  to  use  the  method  of  reflection  just  as  we  did 
for  Laplace's  equation.  For  example,  consider  the  three  dimensional  problem  in 


which  Dirichlet  conditions  are  given  on  a  plane. 
We  reflect  the  source  at  rQ  in  the  plane 
(with  a  sign  change)  and  write 


(t'to) 


-3/2 


-  U(t-t  )  {* 
2a  /v  0 


-|r-ro|2/4a2(t-to) 


Unfortunately,  there  is  no  principle  of  reflection  in  a  sphere  for  the 
diffusion  and  wave  equations. 

3.4  Expansion  in  eigenfunctions 

Let  us  take  a  Fourier  u>- transform  on  t  and  transform  the  equation 

V2G - r  “  -4tt  5(r-r  ) 6 (t-t  ) 

2  it  0  0 

a 

to 

Lot 

V2G  +  ”  G  “  -4tt  6(r-r  )e  0 

a)  2  a)  0 

a 
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Assume  that  ^(r)  are  orthorormal  eigenfunctions,  i.e., 


dt  tp  ,  O' 

n  n 


nn 


2  .2 

corresponding  to  eigenvalues  of  the  problem  ^  +  ft  tp  *  0  plus  whatever  boundary 

conditions  ve  wish  ;o  impose  on  G  (any  weighting  function  can  be  absorbed  by  a 

suitable  definition  of  <l>  ).  The  subscript  n  is  intended  to  be  a  generic 

n 

symbol  which  stands  for  whatever  indices  may  De  present. 

We  substitute  the  expansion 
G 


I  \  *„<*> 


into  the  equation  for  G  : 


“  n  n 

From  the  orthogonality 

a  =  -4tt 
n 

Inverting  the  transform 


o  ,  iwt 

V  a  {-k  +  -~}  o  (r)  ■  -4 r  6(r-r  )e  ° 

2  n  o 

a 


iwt 


(r  )e 


n  o 


2  n 
a 


“  -io(t-t  ) 
o 


G  =*  2ia2  l  (ro)  j  - 

^  w  + 


icA2 

n 


The  integrand  has  a  simple  pole  in  the  lower  half  of  the  to -plane.  We  may  complete 
the  contour  in  the  upper  (lower)  half-plane  for  t-tQ  >  0  .  Thus,  G  vanishes 


for  t-tQ  <0  .We  obtain 


?  2 

.  _  -aV(t-t) 

4ira2  U(t-t  )  l  (r)'l'  (*  )  e 
o  n  no 


In  a  rectangle,  for  example,  with  Dirichlet  conditions  the 
eigenfunctions  are 


<P. 


mn 


2  .  mrx  ,  niry 

-7~r  sin  -  sin  —rd~ 

/ab  a  b 


with 


k2  -  <=42+l^)2 

mn  a  b 
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3-5  Exercises 


1.  Find  the  solution  of  the  one  dimensional  diffusion  equation  for  t  0  which 
satisfies  the  initial  condition 

1|<(X,0)  =  f(x)  _co  <  x  <  05 

2.  Find  the  solution  f or  x  >_  0,  t  >_  0  which  satisfies  iHx,0)  *  f(x)  ip(0,t)  «*  g(t). 

3.  Replace  the  boundary  conditions  of  #2  by<Kx,0)  «  f(x)  ^(0,t)  •  g(t) 

4.  Find  the  solution  of  the  two  dimensional  diffusion  equation  for  t  ^  0  which 
satisfies  the  initial  condition  ij>(x,y,0)  “  f(x,y)  -<»  <  x  <  «  _»<y<»> 

5.  Find  the  solution  of  the  two  dimensional  equation  for  x  >_  0  which  satisfies 

'Kx.y.O)  *  f(x.y)  <KO,y,t)  -  g(y,t) 

6.  Replace  the  boundary  conditions  of  #5  by  tj>(x,y,0)  »  f(x,y)  iji  (0,y,t)  -  g(y»t). 

7.  Find  the  solution  of  the  three  dimensional  diffusion  equation  for  t  >_  0  which 

satisfies  the  initial  condition  4i(x,y,z,0)  »  £(x,y,z)  -a><x<®-«»<y<“ 

_co  <  2  <  “. 

8.  Find  the  solution  of  the  three  dimensional  equation  f or x  >_  0  which  satisfies 

<Kx,y,z,0)  ■  f(x,y,z)  <j>(0,y,z,t)  »  g(y,z,t) 

9.  Replace  the  boundary  conditions  of  #8  by  if*(x,y,z,0)  =  f(x,y,z)  i|>x(0,y,z,t)  « 

g.(y,»z»t)  • 

10.  Find  the  solution  for 

x  ^  0,  y  >_  0,  -co<2<oo> 

which  satisfies 

>Kx,y,z,0)  -  f(x,y,z) 

<H0,y,z,t)  «  g(y,z,t) 

^(x,0,z,t)  «  0 

11.  Find  the  solution  of  the  two  dimensional  diffusion  equation  in  a  rectangle 
of  sides  a  and  b  which  satisfies 

>Kx,y,0)  »  f(x,y)  ^  -  0  on  boundary. 
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Ip(x,b,t)  »  0 


12.  Replace  the  boundary  conditions  of  # 11  by 

<Kx,y,0)  **  4‘(0,y,t)  *  0(a,y,t)  => 

ifrOO.t)  =>  f(x,t) 

13.  Find  the  solution  inside  a  sphere  of  radius  a  which  satisfies 

ip(a,e,<p,t)  =-  f(0,4»,t) 

>Kr>0  «4>>0)  =  0 
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IV.  THE  WAVE  EQUATION 


Green's  function  for  the  wave  equation  satisfies 
2 

f?2 - \  G(r,ro,t,to)  =  -4u  6(r-rQ)6(t-to)  (1) 

c  3 1 

Let  us  try  to  obtain  a  solution  in  the  boundary-free  case  as  we  did  for  the 

diffusion  equation.  We  take  Fourier  k,  ,  k.  ,  k.  transforms  on  x,  y,  z  and  a 

2  ik*r  iut 

Fourier  oj  -  transform:  on  t  obtaining  {-k  +  — ^}G^  =  — 4tt  e  0  e  0 


where,  as  before,  k  =*  (k^.k^jk^) 
function  of  k,  ,  k.,  and  k, . 

1  l.  J 


and  the  subscript  k  indicates  a 


However,  we  find  when  we  try  to  inverts  &ay.,f  the  co-transform, 

2  r  -ico(t-tQ)  ik*r 
,  e  u  o 

dco  — s —  e 


-4ttc" 


2,2  2 

a)  -k  c 


(2) 


f  /C\ 

\ 
/5\  *  \ 

\ 

\ 

- vs - 

-KC- 

+  KC*  J 

The  integral  does  not  exist  since  there  are  poles  at  co  =  +kc  on  the  path 
of  integration  (the  real  axis  of  the  o>  -  plane).  At  this  point  many  writers 

T 

circumvent  the  difficulty  by  saying  that 

the  path  of  integration  must  detour 

around  the  singularities  in  such  a 

way  as  to  satisfy  the  physical  principle  of  causality.  If  the  c  -  function  in 

eq.  (1)  represents  a  pulse  at  r  ,  t  than  we  must  have  G  0  t  <  t 

0  0  o 

since  the  effect  cannot  precede  the  cause.  Therefore,  the  path  is  chosen  as 
shown  in  the  diagram  and  completed  by  a  large  semi-circle  in  the  upper  (lower) 
ha  If -plane  for  t  -  t  >  0  .  The  Green's  function  vanishes  for  t  <•  t0  since  in 

this  case  the  path  encloses  no  singularities. 

There  are  two  criticisms  we  can  make  of  this  approach.  First,  it  is 
not  clear  by  what  authority  one  is  allowed  to  push  the  path  of  integration  off 
the  real  axis.  Secondly,  since  the  problem  is  purely  mathematical  (finding  a 
solution  of  a  differential  equation),  the  resolution  of  the  difficulty  should 
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cal  Ic  is  not  necessary  to  introduce  a  physical  principle, 
ns  leading  to  do  ac  since,  as  we  shall  see,  this  really  has 
causality.  While  it  is  often  helpful  to  think  of  the 
i rt  a  differential  equation  as  a  source  (the  "cause"  of  the 
this  is  not  at  all  a  necessary  interpretation.* 

a..  easily  think  of  the  6  -  function  in  eo.  (l)  as  a 
-nr,:  which  the  fields  are  building  up  throughout 
_  *.  after  which  the  fields  vanish. 


reen'r  function 

-ritr  transforms  has  failed  because  of  a  singularity  on  the 
integral.  Let  us,  therefore,  try  a  Laplace  transform 
:’f  ifilty  cannot  possibly  arise  since  the  inversion  contour 
rt  .f  all  singularities  in  the  complex  plane. 

* is ;  the  inversion  integral  of  the  Laplace  transform 
i  -re.:  the  function  only  for  positive  values  of  its  argument, 
g  zcr  i  for  the  range  (  -  °,0  ) 

f(x)  x  >  0 

X  {f  (X)  }  =  { 

0  x  <  0 


r.  Is  more  reasonable  for  Poisson’s  equation,  for  it  can  be 
'mogeneous  equation  (Laplace’s  equation)  is  satisfied 
e  only  solution  is  the  trivial  one  (a  constant).  There  are 
.•■ources".  But  this  is  not  so  for  the  wave  equation.  There 
iutior.s,  e.g.,  plane  waves,  even  when  the  homogeneous 
led  everywhere  -  effects  without  causes  I  This  is  an  indication 
' re jry  described  by  hyperbolic  equations  the  fields,  e.g.  a 
r  all  of  the  properties  of,  and  are  largely  independent  of, 


Therefore,  since  G  nay  be  taken  to  be  a  function  of  t  -  t0,  we  must  assume  that 
G  (t  -  tQ)  vanishes  for  negative  values  of  its  argument .  This  additional 
condition  is,  of  course,  Just  the  causality  condition  again: 


t  <  t 


(3) 


As  we  shall  see  in  the  next  section,  the  Green's  function  so- defined  is 
waitable  for  a  particular  type  of  initial  value  problem,  viz.,  one  in  which  the 
initial  conditions  are  to  be  projected  into  the  future.  Other  Green's 
functions  can  be  defined  for  other  problems,  e.g. ,  projecting  the  initial 
conditions  into  the  past,  in  which  case  the  causality  condition  is  not 
appropriate . 

Since  G  is  a  function  of  t  -  t0  which  vanishes  for  negative  values  of  this 

argument,  we  may  assume  t  >  tQ  >_  0  .  Taking-  Fourier  ky  k2,  k3 

transforms  on  x,  y,  z  and  a  Laplace  s  -  transform  on  t,  we  obtain  from  eq.  (l) 

*► 

^  0  ik*r  -st 

{-kz-g  /c  }  G.  «  -4ir  e  °  e  0 
ks 

since  G  =  G+  =  0,  t  =  0,  t0  >  0  by  condition  (3). 

Inverting  the  s  -  transform. 


Jk 


.  2  ik*r 

4  TIC  O 

2*i  8 


ds 


J 

<1 


s(t-t  ) 
o 

e  _ 

2,22 

u  +  k  c 


The  integrand  has  simple  poles  at 
(u  =  +  ike  .  By  evaluating  the 
residues  at  these  poles  ve  obtain: 


sin  kc(t-t  )  ik*r 
/  oo 

4rc  - : -  e 


t  <  c 


Now  inverting  the  k. ,  3c„,  k,  transforms  we  have,  for  t  >  t  • 

i  2  o  J 
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G  ■= 


dk^  dk2  dk^ 


^  !J 


8 in  kc(t-t  )  -  ilc-{  r- r  ) 

- c - 2-  « 


The  integral  is  evaluated  by  means  of  the  following  trick.  We  introduce 
spherical  coordinates  in  the  k^,  k^,  space  so  that 


dk^  dk,  dk^  a  k  sin0  d6  d4>  dk 


and  we  orient  the  coordinate  system 

so  that  6  is  measured  from  the 

direction  of  r-r  .  The 

o 

cp  integration  can  be  performed 
immediately,  giving 


4ttc 

(2*r 


kdk  sinkc  (t-t  )  d9  sin8  e 
o  ' 


-ikjr-ro JcosS 


0  0 
The  6  integration  is  also  readily  performed: 

-ik| r-r  jcos8  2  sink  jr-r 


d6  sinG  e 


J 

0 


k  | r-r  | 
1  o' 


We  express  the  sines  as  exponentrals ,  and  since  the  integrand  is  even  we  may 
write 


-c/4tt 


-*■  -h 

r-r 


ikc(t-c  )  -ikc(t-t  ) 
dk{e  °~e  °} 


o  5  -® 


ikjr-r  |  -ikjr-r  | 

,  1  o'  1  o' , 

te  -  e  } 


Recalling  the  vellknown  representation  of  the  delta  function 


«(p) 


J. 

2tt 


a1**  dx 


ve  have 


G  * 


-c/2 


r-r 


+  c(t-tQ)] 


+  6[-|r-ro|  -  c(t-to)] 

-  6[|r-?o|  -  c(t-t0)] 

-  6 [- 1 r-?o j+  c(t-to)]  } 

Since  we  have  assumed  t  >  t  >  0  the  arguments  of  the  first  two  delta  functions 
can  never  be  zero  and,  therefore,  these  functions  vanish.  The  other  delta 
functions  are  identical  because  the  Function  is  even.  Making  use  of  the 
identity 


<5  (ax) 


<$(x) 


we  have 


G 


r-r 


(3a) 


We  have  obtained  this  expression  assuming  t  >  tQ  >  0  #  However,  we  now  notice 

that  the  expression  automatically  gives  zero  for  t  <  tQ  an<*»  therefore,  is  correct 
for  all  t,  tQ. 

If  we  think  of  G  as  the  response  to  a  pulse  at  r  ,  t  ,  then  we  see  that 
this  response  is  itself  a  (spherical)  pulse  expanding  with  a  velocity  c  and 
damped  by  the  reciprocal  of  its  radius. 

In  the  boundary- free  case  with  no  initial  conditions  (i.e.,  assuming  the 
sources  vanish  sufficiently  far  in  the  past)  a  solution  of  the  inhomogeneous 
wave  equation 

2 

{V2 - -  =  ~4tt  f(r,t) 

c2  at4 
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May  be  written 


iKr.t) 


dt 


dt  f(r  ,t  )  G(r,r  ,t,t  ) 
o  o’  o  ’  o’  o 


and  using  the  sifting  property  of  the  delta  function  to  perform  the  t  integration 
we  have 


dT 


"  IK  l 


This  solution  is  sometimes  called  the  retarded  solution.  It  is  an  integration 
over  the  sources  evaluated  not  at  t  but  at  a  certain  time  before  t,  the  retarded 
time  t  -  jr-ro|/c  .  The  dependence  on  the  retarded  time  reflects 

■4* 

the  finite  velocity  of  propagation.  Any  change  in  the  source  at  vq  is  not 

felt  at  r  until  a  time  jr-r  |/c  later. 

o 1 

Let  us  again  emphasize  that  the  above  is  an  interpretation  of  the  purely 
mathematical  problem:  find  a  solution  which  vanishes  as  t  -«  of  the 

wave  equation  with  an  inhomogeneous  term  which  is  non- zero  in  a  finite  region 
of  space-time.  The  interpretation  consists  of  identifying  the  inhomogeneous 
term  as  the  source,  or  cause,  of  the  solution.  We  can  also  find  a  solution 
which  vanishes  as  t  -*■  «  by  using  a  different  Green's  function,  the  so-called 
advanced  Green's  function  (see  §4.3). 

4.2  The  initial  value  problem 

We  now  assume  that  we  seek  a  solution  for  t  >  0  which  satisfies  initial 
conditions  at  t  =  0  as  well  as  boundary  conditions  on  whatever  surfaces  may  be 
specified.  The  wave  equation  is  symmetric  with  respect  to  past  and  future,  but 
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the  statement  of  the  initial  value  problem  introduces  asymmetry:  we  wish  to 
continue  the  initial  conditions  into  the  future  rather  than  into  the  past.  The 
asymmetry  is,  for  thi3  problem,  appropriately  introduced  into  the  Green's 
function  by  imposing  the  causality  condition  (3) . 

We  first  establish  the  reciprocity  relation 

G(r,r  ,t,t  )  *  G(r  ,?,-t  ,-t) 

OO  0  0 


The  relation  between  t  and  t  would  hold  for  any  function  of  t  -  t  ,  while  the 

o  J  o 

h*'  -y 

relation  between  r  and  rQ  is  stronger,  expressing  the  equivalence  of  source  and 
field  points.  We  first  write 


- "  -4"  <5(r-r  )<5(t-t  ) 

c  3t2  0  °  o  o 

{V“ - j  -22’-  G(r,r1>-t,-t1)  =  -4tt  dCr-r^dCt-t.,) 

c  3 1 

We  multiply  the  first  equation  by  GCr.r^-t.-t^) ,  the  second  by  G(r,rQ,t,to) , 
subtract  and  integrate  over  t  and  whatever  volume  is  of  interest. 


CO 


—00 


dt 


dr  {G(r,r1,-t,-t.)V2  G(r,r  ,t,t  ) 
11  o  o 

G(r,ro,t,tQ)V2  G(r,r1,-t,-t1)} 


dt 


dt  ~  {G(v,vv- t,-tl)  g|c(?,?o,t,to) 


-  G(r,rQ,t,to)  GCr.r^-t.-tjH 


"  4"  VW  "  G^o’^l,-to,-tl^ 
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3y  Green's  theorem  the  first  integral  is 

o© 

J  dt  \  ds  {GfJ.V-t.-V  ^  G(?,ro.t,to) 


"  G(r*fo‘t‘to)^G(r‘rl>-t’-tl)} 


We  now  impose  the  condition  that  G  satisfies  homogeneous  boundary  conditions  of 
the  type  given  on  the  surface.  For  example,  if  Dirichiet  conditions  are  given 
G  must  vanish  on  the  surface;  if  Neumann  conditions  are  given  c> G/^n  must 
vanish  on  the  surface.  In  the  latter  case  there  is  no  arbitrary  additive 
constant  as  there  was  for  Laplace's  equation  because  any  such  constant  vr’ild 
be  fixed  by  the  initial  conditions.  Having  imposed  this  additional  condition 
the  surface  integral  vanishes.  The  second  integral 


f  dt 


[G(rtr  ,-t  ,-t, )  -r—  G(r,r  ,t,t  ) 

i  x  at  O  O 


-  G(rtSo,t,to) 


(*0 


vanishes  because  of  the  causality  condition  and,  therefore,  the  theorem  is 
proved. 

To  solve  the  initial  value  problem  we  write 

9 


-7 -^o-}  ,t  )  = 

c2  at2  0  0 
% 

-r  — r)  G(r,r  ,t,t  ) 
c2  at2  0  0 


— Att  f(r  ,t  ) 
o  o 

=  —Att  6(r-r  )6(t-t  ) 
o  o 


o 

The  second  equation  follows  from  eq.  (l)  by  the  reciprocity  relation.  We 
multiply  the  first  equation  by  G,  the  second  by  'P  (  r0»t0  ),  subtract  and 
integrate  over  the  volume  of  interest  and  over  0  t  <  “  .  Now, 


CD 


-  G 


ii-i" 

at  Jo 

o 
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which  vanishes  at  the  upper  limit  because  of  the  causality  condition.  Using 


Green's  theorem 

f 

Kt.t)  -  j  dto  j  dTo  f(ro,to)  G(rfr0,t,t0) 
0 

c+ 


4tr 


f 


dt  i  ds  {G  -  *  I5-} 
o  j  o  3n  3n 

J  oo 


4itc 


,  ,r  ,  3G  , 

dTo  [G  n  '  ♦  ?r!t  -o 

O  0  0 


(5) 


Because  of  the  causality  condition  the  integrals  over  t  have  reduced  to  integrals 
over  [o,  t  *■]  where,  as  before,  t  +  indicates  t  +  e  in  the  limit  e  0+ 

Eq.  (5)  is  a  solution  of  the  boundary  value  problem.  The  boundary  conditions 
enter  through  the  surface  integral  in  the  second  term.  We  recall  that  we  have 
assumed  that  G  staisfies  homogeneous  boundary  conditions  of  the  type  imposed  on 
.  If  Dirichlet  conditions  are  given  then  G  vanishes  on  the  surface  and  the 
integrand  is  i  3G/3nQ 

If  Neumann  conditions  are  given  the  integrand  is  G  3t|>/3no  .  If  mixed 

conditions  specifying  +  a  3^/3n  on  the  surface  are  imposed  then  G 

satisfies 

3G 

G  +  a  ~  »  0 
3n 


on  the  surface  and  the  integrand  may  be  written 


G 


2£_ 

3n 

0 


—  G  +  } 

a  3n 

0 
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The  effect  of  the  initial  (Caichy)  conditions  which  specify  $  and  its  time 

derivative  at  t  =  0  is  contained  in  the  third  integral  in  eq.  (5). 

4.3  Other  Green's  functions 

The  Green’s  function  of  the  last  two  sections  is  sometimes  called  the 

retarded  C-reen's  function,  G^.  If,  given  the  initial  conditions,  we  wish  to 

find  the  solution  at  some  time  in  the  past  G  .  is  no  longer  appropriate. 
Instead  we  must  use  the  advanced  Green's  function,  >  which  is  obtained 

by  imposing  the  condition 


G  -  0  t  >  tQ  (6) 

which  we  might  call  the  "anti- causality  condition".  We  see  that  auxiliary 
conditions  such  as  (3)  and  (6)  are  only  devices  used  to  obtain  solutions  of 
certain  types  of  problems,  and  need  not  be  given  a  causal  interpretation. 
Indeed,  in  the  case  of  eq.  (6)  such  an  interpretation  would  say  that  the  effect 
precedes  the  cause! 

In  the  boundary- free  case  G  ,  can  be  calculated  by  the  methods  of  §  4.1 

fluV 

giving 


r-r 


adv. 


r-r 


■  0  ■ 

The  initial  value  problem  is  solved  as  in  §  4.2,  integrating  t0  over  [-  •*  ,  0]. 

Finally,  we  note  that  G  ,  .  satisfies  the  same  reciprocity  relation  as  g  •'  for 

adv .  ret . 

condition  (6)  also  causes  the  integral  (4)  to  vanish. 

A  Green's  function  which  is  symmetric  in  the  time  can  also  be  constructed. 
We  have  noted  the  difficulty  associated  with  the  integral  (2)  due  to  the  simple 
poles  on  the  path  of  integration.  By  using  the  theory  of  distributions  it 


a 

? 


a 
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can  be  shown  that  it  is  correct  to  take  the  principal  value  of  this  integral 
(Light hi 11,  1958,  p.3l).  It  is  readily  shown  that  the  Green’s  function  obtained 
in  this  way,  G  say,  is 


G 


+  G 


adv. 


} 


This  Green's  function  is  useful  in  problems  which  are  symmetric  in  the  time,  for 
example,  in  Dirac's  theory  of  the  electron  and  in  certain  problems  in  quantum 
field  theory.  See  also  §  5*3* 


4.4  Method  of  reflection 


We  turn  now  to  the  problem  of  constructing  Green's  functions  which  satisfy 
homogeneous  boundary  conditions  on  surfaces.  For  the  wave  equation  the  method 
of  reflection  is  unfortunately  limited  to  conditions  on  planes.  For  example, 
suppose  we  seek  a  solution  in  the  half- space  x  >  0  with  Dirichlet  (Neumann) 
conditions  on  the  plane  x  »  0.  We  must  construct  a  Green's  function  which 
satisfies  G  =  0  (  3G/3n  =  0  )  on  x  «  0  . 

This  can  be  accomplished  by  placing  ^  , _ ^ 

image  at  v±  =  (  -  3^,  jf0,  zQ)  which 


an 


is  the  reflection  of  rQ  in  the  plane. 


O 


G 


I-*  "*■ 

r-r 


I-*-  -* 

r-r 


SU-to  - 


-)  + 


r-r. 


$  (t-t 


The  upper  sign  is  used  for  Dirichlet  conditions,  the  lower  for  Neumann. 


If  we  think  of  the  first  term  as  a  spherical  wave  spreading  out  from  the 


point  ,  then  we  may  think  of  the  second  term  as  a  wave  reflected  from  the 
plane. 
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k .  [>  Method  of  eigenfunction  expansions 

This  method  is  more  general  and  more  complicated  than  the  method  of  reflection. 


We  operate  on  eq.  (l)  with  a  Laplace  s- transform  on  t  obtaining 

?  9  -►  -y  -8t 

(V  +  k  )G  »  -4tt  6(r~r  )e  0  t  >  0 

so  o 

where  k  =  i  s/c.  Let  (r)  be  ortho-nomal  eigenfunctions  of  the 

n 

homogeneous  equation 

(92  +  k2)^  -  0 

n  n 

(the  Helmholtz  equation)  which  satisfy  the  boundary  conditions  imposed  on  G. 
For  the  eigenvalues  k  and  the  eigenfunctions  <L  the  subscript  n  is  a  generic 
symbol  which  stands  for  whatever  indices  may  be  present. 

For  example,  for  the  interior  of  a  unit  sphere  the  (unnormalized) 
eigenfunctions  which  vanish  on  the  surface  are 

Wr'M>  '  Jtj.  i  <V)P»tcos9)  elm* 

2 

where  the  eigenvalues  are  the  roots  of 


which  can  be  shown  to  be  real.  Here  1  and  m  are  integers  and  the  subscript  n 
is  used  to  number  the  infinitely  many  roots  of  the  1^  equation. 

We  substitute  the  expansion 


(7) 


into  eq,  (7)  obtaining 


£  a  {-k2  +  k2}i^  (r)  =  -4*  6(r~r  )e 


-st 


n  n 


and  since  by  our  assumption  of  ortho-normalized  eigenfunctions 


dr  ’4  (r)0  , (r)  =  6  , 


n  n 


nn 


we  obtain 


a  =  -4t 


*  -*•  _sto 
>P  (r  )e 
n  o _ 

2  2 
k  -k 

n 


=  4irc 


2  Wc 


"St 


0  2  2 

s  +  k  c 
n 


Inverting  the  s -transform 


„  .  1  G  eSt  ds 

2iri  s 


(r)ip  (r  )  . 
2tt  i  Yn  n  o  J 


s (*■— t  ) 

0 

e _ 

2.22 
s  +  k  c 
n 


ds 


where  the  contour  c  is  to  the  right  of  all  singularities  (in  our  example  of 
the  sphere,  these  are  on  the  'imaginary/  axis  of  the  s  -plane). 

Evaluating  the  residues, 


lk»c(t-to)  -lknc(t-to) 

G  -  4,0  I  *n(b*n(r  >  2 - 2TTT - 

n 


sin  k  c(t-t  ) 


4itc  £ 


(r) 

non 


4.6  Exercises 

1.  Obtain  the  boundary- free  retarded  and  advanced  Green's  function:: 
for  the  one  dimensional  wave  equation: 
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1  -  U ( : x-x  !  +  c( t-t  ) ) } 
0  o 


x  >  0 


x  <■  0 


rtair.  Alembert's  solution  of  the  initial 


.  <  x-ct  ,0)  +  »  (x+ct  ,0)  +—  1  i,  U,Q)dO 

c  ;  t 


va! ;.{  for  both  positive  and  negative  t  Hint: 
i  u.l*.  function  is  a  delta  function,  dL'/dx  - 


r.e  4if.eri3ior.al  wave  equation  for  x  >  0,  t  >  0 


i\xi  ,r(x,0)  =  g(x)  t*  ( 0 ,  t )  »  h(c) 


•  .or..'  of  “3  by 


x;  ,  (x,0)  »  g(:<)  tx(0,t)  »  h(t) 


,'ie:  :  O' ir.d&ry- free  Green's  function  for  the  two 
•-q.f-tioi.  is  given  by 


r-rol/c} 


a- race  transforms  lead  to  a  difficult 


:  "ir.t  soui'ce  in  two  dimensions  is 


urce  in  three  dimensions.  Therefore, 


iafttuUnvixtaM  ajaSinikliisui 


G  can  be  obtained  from  the  *  ~ee  dimensional  Green’s  function 
by  integration  over  R0- 

6.  Solve  the  initial  value  nro'ciec  in  two  dimensions  assuming 


that  no  boundaries  are  oreser.t. 


7.  Find  a  solution  of  the  three  d isiets.  ;nsl  wave  eouatior.  .or 


z  >0  which  satisfies 


u(x,y,G,t)  “  f (x,y , t) 


-=£  <.  c  <.  X 


Eigenfunctions  such  as  those  giver,  in  S  4 . £  for  the  Interior  of  a 
sphere  are  not  always  available.  For  example,  for  the  exterior 


of  a  stthere  ve  reauire  functions  which  give  or.lv  outgo  in 


at  infinity.  This  is  the  so-called  radiation  condition;  it  is 
equivalent  to  the  causality  condition.  Therefore,  instead  of 

the  Bessel  function  we  would  have  to  use  H^.(kr  ) 

(!)  i+  1 
But  H'~'  r.as  r.o  reaj.  zeros  ar.c,  therefore,  cannot  satisfy  the 

boundary  condition  on  tno  oonsnS’ . 

Vhen  an  expansion  of  the  fundamental  solution  (such  as  that 

used  ir.  problem  2.12 )  is  known,  the  following  approach  may  be 


<1 


used.  Let  Gv  =  Gs  exp  stQ .  Show  that  in  the  boundary- free  case 


G.  -  ^eikR 
k  R 


where  R  =  |r-r 


a 


so  that  generally  G>  behaves  as  (exp  i  k  ?. )  /?.  near  the 
singularity  at  rQ  . 

Therefore,  write 


+  F(r) 


m 


-  kh  - 


where  7  (  r  )  is  regular  in  the  region  of  interest.  By  expanding 


”  in  appropriate  eigenfunctions  and  making  use  of  the  expansion 
(Morse  and  Ingard  .  p.  352) 


ikR 


=  ik  :  (2£+l) 


<H 

[raj 

ILL 

(£+ 

ft] 

!)! 

cor  struct  a  1.  f\. r  the  exterior  of  a  snhere  of  radius  a  which 


vanishes  or.  the  surface.  In  the  above  expansion  j  and  are 

spherical  Bessel  and  Hansel  functions,  r>  (  r<  )  is  the  larger 

(smaller)  of  r,r  ,  ar.d  Y  .  (  0,b  )  -  7,111  /  - — Q 

o  m*.  ' 


(  cosG  )  e 


V.  ACOUSTIC  RADIATION  AND  SCATTERING 


Our  purpose  in  this  chapter  is  not  to  give  a  general  introduction  to 
acoustics,  but  simply  to  give  some  examples  of  the  use  of  the  Green's  function 
for  the  scalar  wave  equation.  Readers  seeking  a  general  background  in  acoustics 
or  examples  with  more  complicated  physics  than  those  presented"  here  are  refered 
to  Morse  and  Ingard  (1968). 

5*1  Basic  Equations  of  Acoustics 

For  simplicity  we  assume  a  homogeneous,  inviscid,  compressible  fluid. 

The  conservation  of  mass  is  expressed  by  the  continuity  equation 

dp  , 

—  -r  dlv  pu  -  q 
0  Z 

where  0,  p,  and  q  are  the  density,  velocity,  and  source  strength  (rate  at  which 
fluid  is  "created"  (introduced  into  the  flow  field  from  outside)  per  unit  volume. 
The  force  equation  (Newton's  second  law)  for  fluids  is  usually  called  Euler's 
equation: 

p  +  pp*?U  u  —  Vp  +  F 


where  p  is  the  pressure  and  F  the  external  force  per  unit  volume.  The  field 
variables  (unknowns )  are  p,  p,  and  ,1 .  Therefore,  another  relation  is  needed 
to  complete  the  above  equations.  For  this  we  assume  adiabatic  motion  (which  is 
quite  realistic  for  sound)  so  that 


p  =  const,  p 

where  y  =  c  /c  ,  the  ratio  of  specific  heats  at  constant  -pressure  and 
p  v 

constant  volume. 

The  basic  equations  of  acoustics  are  obtained  by  linearization. 


We  write 


p  =  0  +  p' 

o 


p  =  p  +  p' 

■  o 


F  =  F  +  F' 
o 


u  =  u  +  m' 
o 


+  q' 


where  the  quantities  with  subscripts  are  mean  or  equilibrium  values,  which 
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satisfy  the  equations  separately,  and  the  primed  quantities  are  fluctuations 


assumed  small  with  respect  to  mean  values.  We  assume  that  there  is  no  mean  flow 
-*■ 

so  that  mq  =  0  and  substitute  these  relations  into  the  above  equations 
neglecting  quadratic  terms  in  the  primed  variables.  Prom  the  adiabatic  law 
we  oo>.ain 


P’ 


2 

=  c  p 


f 


i  /  2 

Where  c  =>  (yp  / p  ) "  is  the  velocity  of  sound.  Using  this  relation  to 

oo 

eliminate  p',  we  obtain  by  linearization  of  the  first  two  equations  the 
basic  field  equations  of  acoustics: 


1  5p  .  . . 

1  at  +  ?o  dlv  u  a  ’ 


%  Tt  + 


(1) 

(2) 


where  for  convenience  the  primes  have  been  dropped  from  perturbation  quantities. 

We  now  note  two  consequences  of  these  equations.  By  multiplying  the 

first  by  p,  the  second  by  and  adding  we  obtain  the  conservation  of  energy 
3w  -*• 

■r-  +  div  s  =  pq/p  +  F*y 
o  t  O 

where  w  is  the  energy  density 

1,2,  2  ,  2. 

w  =  -  (p  /pc  +  p  u  ) 

£.  O  O 

and  s  (the  Poynting  vector)  is  the  energy  flux  density  (energy/ cm^  sec  ) 


The  terms  on 

pressure: 


S  =  py"’  (3) 

the  right  hand  side  represent  dissipation. 

eliminating  y  between  eqns.  (l-2)  we  obtain  a  wave  equation  for  the 

.  r.J-la 


5-2  Multipole  Fields 


We  will  not  be  concerned  with  the  mechanism  of  radiation  in  this  section, 
but  will  simply  assume  that  there  is  a  finite  source  region,  i.e.,  a  region  in 
which  energy  is  somehow  introduced  into  the  acoustic  field.  The  mechanism,  for 
example,  could  be  a  vibrating  surface  or  turbulence.  We  further  assume  a  har¬ 
monic  time  dependence.  A  typical  time  dependent  quantity  is  then  written 

p(r,t)  =  p^(r)e 

This  assumption  is,  of  course,  equivalent  to  tds ing  a  Fourier  transform. 

We  write  the  source  term  in  the  wave  equation  as  a  function 
f (r,£)  «  f  (r)e 

Then  the  wave  equation  (4)  for  p  becomes  the  Helmholtz  equation  for  p  : 

2  2 

V  p  +  k  p  »  -4*  f  k  «  w/c 

Solving  the  wave  equation  by  means  of  the  boundary- free  Green’s  function  (4.3a) 
gives 


-iwt  j 
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J 

0  0 

ik 1 r-r  1 

w 
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f  (r  ) 
w  0 
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dr  G  (r,r  )  f  (r  ) 
o  w  o  wo 


where  is  the  boundary- free  Green’s  function  for  the  Helmholtz  equation. 


G  (r,r  ) 
w  o 


exp  ikjr-r^ 
r-r  | 


(5) 


We  assume  that  the  origin  of  our  coordinate  system  is  somewhere  with¬ 
in  the  souce  region  and  that  we  wish  to  calculate  the  field  far  from  this  region 
so  that  £  is  large  compared  with  the  values  of  rQ  within  its  range  of 


integration.  We  expand 


i-* 

r-r 


in  a  Taylor  series  in  about  the  origin 


l"v  "v 

l  r-r. 


r’r 


=  r  - 
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Then  the  expansion 


ik-r-r 

o 


~  e 


ikr  /t  ,  i 

i  exp  ikr^  •  {—  +  . . . } 

ikr 


o  r 
■  r 


11  +  ik 


,} 


is  valid  provided  we  make  the  further  assumption  kr  <<  1  .  This  means  that 

*  o 

we  assume  the  dimensions  of  the  source  small  compared  with  the  wavelength  * 


T 

(k  =  2^ /> )  .  with  this  assumption  we  may  expand  in  a  Taylor  series  in  r0  , 

I 


eikr  d 

T~  +  4 
1 


ik|r-rQ! 


I I  r  a0  Xoi  + 
oi  I r-r  I  o 


Substituting  into  the  integral 


we  obtain 


dr  f  (r  )G 
I  O  0)  0  to 


ikr 


ikr 


“w  r  r 

where  we  might  call  S  the  monopole  strength 


+  .  V  . 


S  = 


dr  f  (r  ) 
O  u!  O 


since,  as  we  see  from  the  Green's  function,  the  first  term  is  the  field  of  a 
point  source  at  the  origin. 

In  writing  the  second  term  we  have  used  the  fact  that  3  ”^/,3xi 

when  these  operate  on  functions  of  lr-r0i  only.  The  second  term  is  a  dipole 
field  with  dipole  moment 


D.  = 

l 


dt  x  f  (r  ) 
o  oi  w  o 


Higher  order  terms  in  the  expansion  represent  quadrupoles,  octupoles 
etc.  These  terms  reflect  the  symmetries  of  the  source.  For  example,  a  sphere 
vibrating  purely  radially  has  only  a  monopole  term,  which  is  the  only  spherically 
symmetric  multipole.  The  dipole  term  reflects  symmetry  about  an  axis  and  the 
higher  order  multipoles  are  correspondingly  more  complicated. 
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5.3  Radiation  Calculations 

V/e  will  discuss  two  methods  of  calculating  the  energy  radiated  by  sources. 
The  first  and  best  known  of  these  we  call  the  far  field  method.  It  consists  of 
integrating  the  Poynting  vector  s  (3)  over  a  sphere  whose  radius  is  arbitrarily 
large.  Thus,  it  represents  only  energy  which  is  ultimately  lost  by  the  system, 

i.e.,  radiated  away.  Since  the  element  of  surface  area  in  spherical  coordinates 

2  *+ 
is  r  sinG  d0  d<f>  where  r  is  the  radius  of  the  sphere,  and  since  s 

is  quadratic  in  the  fields,  only  terms  of  order  1/r  will  contribute  to  this  integral 

in  the  limit  r  -*■  »  .  These  terms  are  called  the  radiation  fields. 

For  example,  the  field  of  a  dipole  along  the  z  -  direction  is 

i(kr-cct) 


d_ 

dz 


z  z  .  i(kr-mt) 

=  (ik  — - jJe 

r  r 


The  first  term  gives  the  radiation  field  since  it  is  the  only  term  of  order  1/r. 


rad. 


ik  cos0  i(kr-ut) 

-  e 

r 


The  velocity  u  in  the  radiation  field  is  obtained  from  eq.  (2) 


1  "  po  “rad.  *  prad.  (ikz/r) 


and  the  time-averaged  Poynting  vector  is 

s 


rad. 


<Re  p  ,  Re  y  ,  > 
~  rad.  rad. 

1  k  cos  ^6 


2  p  c  2 
o  r 


in  the  radial  direction.  This  shows  the  directional  properties  of  the  radiation. t 
The  total  power  may  be  obtained  by  integration. 

In  addition  to  the  radiant  energy  in  the  field  there  is  also  reactive 
energy,  energy  which  fluctuates  between  the  source  and  field.  The  second 
method,  the  near  field  method,  consists  of  integrating  the  energy  flux  density 
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over  the  source,  nils  has  the  advantage  of  giving  the  reactive  as  well  as  the 


radiant  energy,  and,  therefore,  the  mechanical  impedance. 

For  example,  assume  the  normal  velocity  y(r)e  is  specified  on  a 

vibrating  surface,  Py  eq.  (2 )  this  is  equivalent  to  giving  the  normal  derivative 
of  the  pressure  (we  are  still  assuming  harmonic  time  dependence).  Assuming  the 
(retarded)  Green's  function  is  known  we  may  substitute 


JL 

Ait 


G  ds  dt 
3n  ret.  o  o 


iwp^ 

An 


p(r  )e.  iu)t  G  dt  ds 
o  ret.  o  o 


into  the  expression  for  the  total  power 


P  *»  j  p  y  ds 


(6) 


We  may  now  write 


=  G  +  H 


ret. 

ret. 


+  G 


adv. 


-  G 


adv. 


} 

} 


where  G  is  the  time  symmetric  Green's  function  of  §4.3*  H  is  not  a  Green's 
function  since  it  is  a  solution  of  the  homogeneous  wave  equation,  as  we  can  see  by 
'subtracting  the  equations  for  Grgt  and  . 

After  substituting  into  eq.  (6)  we  find  that  if  the  sign  of  t  is  changed 
the  part  involving  G  changes  sign  and  the  part  involving  K  is  unchanged  Therefore, 
the  former  represents  reactive  power  and  the  latter  resistive  or  radiant  power. 

This  elegant  method  was  perhaps  first  used  by  Schwinger  (1949)  in  connection 
with  electromagnetic  radiation. 

As  an  example  consider  an  infinite  plane  z  =  0  on  which  y  =  0  except  for  a 
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small  portion  0  which  vibrates.  The  Green’s  functions  are  easily  obtained  by 
the  method  of  reflection  with  an  image  ppint  *  (xo»yo>-zo)  ‘ 

effect  of  the  harmonic  time  dependence  is  to  replace  G  ret  a(^v  .  by 

+ik | r-r  q |  +ik | r-r^ j 

^wret. ,adv. 

On  the  plane  lr-r0l  “  } r ^ 1 
radiative  power 


r-r  r-r. 

o'  1  i 

.  We  have  for  the  time- averaged 


rad. 


u(r)  ds 


top 


~  u(r  ) 
2tt  o 


sin  k|r-r  j 
_ 1  o' 

I  r-r. I 


ds 


and  for  the  reactive  power 


top 


o 

2u 


ds 


cos  k  r-r 


ds  y(r)  y (r  ) 
o  o 


i  **■ 

r-r 


(7) 


The  integrand  in  the  last  equation  is  singular,  but  by  a  well  known  theorem  for 
improper  integrals  (Kellogg,  p.  lb-9)  the  integral  is  convergent. 

Detailed  calculations  are  available  for  the  case  in  which  0  is  circular. 

If,  further,  y  is  constant  (circular  piston)  then  the  coefficients  of  in  the 
expressions  for  Pra(j  .  and  Px  give  the  resistive  and  reactive  parts  of  the 
complex  mechanical  impedance.  For  further  details  see  Bouwkamp  (1946)  and 
Morse  and  Ingard  (p.  38l  ff.). 

If  directional  properties  of  the  rauiation  are  desired  we  can  make  use  of 
the  relation  (  §  4.l) 


r  ikn* (r-r  ) 
e  0  dQ 
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where  d.: 


the  element  of  solid  angle  about  the  unit  vector  n  .  Then  from  eq. 


(7) 


where 


^rad. 


dn  p(0,<J>) 


up  k 
o 


-v  -y 

r  -ikn*r  , 

{  ds  u(r  )e  } 

J  o  o 


is  the  power  radiated  into  dft  .  If  0,<j>  are  spherical  angles  with  the  normal 
to  the  plane  as  the  z  -  direction  and  4>  is  the  azimuthal  angle  in  the  plane, 

then 

n*r  =  r  cos  (n,r  ) 
oo  o 

■-  r  sin0  cos  (d>~ d>  ) 
o  o 


in  the  above  formula. 

5 . U  Mechanisms  of  Radiation 

Let  us  n^w  take  a  closer  look  at  the  source  terms  in  eq.  (k).  We  shall  add 

to  these  a  term  from  the  stress  te  sor  which  is  of  the  order  neglected  in  the 

derivation.  Lighthill  (1952, 1954)  has  shown  that  this  term  may  not  be  of 

second  order,  at  least  for  turbulence  in  air.  Therefore,  we  write 

2 

V2p  +  k^p  =  iuq  +  7-F  -  l  .  r—  U  .U  . 

*u  *u  nu  w  dx^oXj  o  ui  uj 

Solving  by  means  of  the  Green's  function  (5)  we  have 


di  G  (r,r  )  {iuq  +  V  *F  -  £ 

o  u  o  u  o  u 


3x  ,9x  .  PotJui"uj 
oi  oj 
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(3) 


■} 


where  3/3x  f  * 

Let  us  look  at  the  terms  in  this  integral  separately.  If  we  expand  the 
first  integral,  the  one  containing  ,  as  in  §  5.2  the  first  term  in  the 
expansion  will  be  a  monopole  field.  For  this  reason  (or  3q/ 3t)  is 

called  a  monopole  source.  It  represents  the  rate  at  which  fluid  is  added  or 
withdrawn.  A  vibrating  surface  is  a  source  of  this  type,  for  its  effect  is 
as  if  fluid  were  alternately  injected  and  withdrawn.  An  example  is  tne 
cavitation  behind  a  pi^p^ler.  A  collapsing  babble  is  cushioned  by  vapor  inside 
it  causing  a  bouncing  effect.  Of  course,  «...  "hould  keep  in  mind  that  the  field  is 
not  a  pure  monopole  field,  but  rather  the  monopole  term  is  the'  -dominant  term  in  a 
multipole  'mansion. 

The  second  integral  may  be  transformed  through  the  identity 


—y  ~y 

~y 

dr  7  ‘G  F 

dr  G  V  *F  + 

0  0  CO  (0 

3 

O 

3 

O 

dr  F  *V 
O  (0  o 


G 

CO 


We  have  assumed  that  the  sources  are  non- zero  in  a  finite  region.  The  left  hand 
side  may  be  written  as  a  surface  integral  by  Grauss1  theorem,  but  since  the 
volume  integral  is  over  all  space  the  limit  must  be  taken  as  the  surface  goes  to 
infinity.  As  the  limit  is  taken  the  surface  will  eventually  lie  in  a  region 
where  F  vanishes.  Hence,  the  left  hand  side  is  zero  and  the  contribution  to 

(i) 

the  pressure  is 

dr  F  *V  G 

0  (0  0  CO 

T1  a  leading  term  of  a  multipole  expansion  of  this  integral  will  contain 

-V 

derivatives  of  ;  it  is  a  dipole  term  and,  therefore, 


is  called  a 


dipole  source.  An  example  is  the  sound  field  due  to  the  force  exerted  or  &  fluid 
by  a  propeller  blade. 

-The  third  integral  in  eqs  (8)  can  be  transformed  two  integrations  by 
parts  giving 

32 

dr  I  p  p  .11  .  ~ — s — “  G 

J  o,i  ..l.j  »*0l»*0j  « 

The  leading  term  in  a  multipole  expansion  contains  second  derivatives  of  Gu; 
consequently,  it  is  a  quadrupole  term.  Thus,  for  example,  a  turbulent  region  is  a 
source  of  quadrupole  radiation. 

5*5  Scattering  and  Diffraction 

Scattering  and  diffraction  refer  to  the  same  phenomenon,  the  interaction  of  a 
wave  with  an  obstacle.  The  word  diffraction  is  usually  used  when  the  dimensions 
of  the  obstacle  are  large  compared  with  the  wavelength,  as,  for  example,  when  sound 
passes  through  an  -^erture  in  an  infinite  screen.  The  obstacle  need  not  be  a  solid 
body.  It  may,  for  example,  be  a  region  of  turbulence  or  inhomogeneity. 

Green’s  function  is  used  in  these  problems  to  transform  the  boundary  value 
problem  into  an  integral  equation.  Let  us  consider  the  case  of  a  solid  body. 

For  harmonic  time  dependance  and  with  no  sources  present  a  solution  is  given  by  the 
surface  integral 

-  r  3p  3G 

j  d,0  G»  - 

°  °  (9) 

This  is  not  a  solution  of  the  boundary  value  problem  since  it  is  not  proper  to 
specify  both  Py  and  its  normal  derivative  on  the  surface.  However,  instead  of 
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new  imposing  s  boundary  condition  on  0  as  ve  |»v6  in  previous  proble&a,  ye  my 
instead  regard  eq»  (9)  as  an  integral  equation  for  Pu .  The  solution  is  not 
unique  since  we  have  not  yet  imposed  any  boundary  conditions.  We  may  add  to  it 
any  solution  of  the  homogeneous  wave  equation.  Therefore,  we  write 

p  -  Pj.  +  Jr  )  ^3^  -  Poa-T  \ 

where  is  a  given  incident  wave.  The  integral  represents  the  scattered  wave  if 
we  use  the  retarded  Green’s  function 


* 

which  gives  outgoing  waves.  The  boundary  condition  may  now  be  imposed.  For  a 
rigid  body  we  would  set  the  normal  derivative  of  the  pressure  equal  to  zero  on 
the  boundary.  For  a  compliant  body  a  linear  combination  of  and  its  normal 
derivative  would  be  specified. 

For  simple  geometries  the  integral  equation  can  often  he  solved  by  expanding 

0  ,  p  ,  and  G  in  appropriate  orthogonal  functions,  fixing  undetermined  coefficients 

*ui  ri  w 

by  evaluating  on  the  surface.  Such  problems  can  be  solved  by  direct  methods, 
so  that  nothing  is  really  gained  by  the  transformation  to  an  integral  equation. 

When  an  exact  solution  is  not  possible,  however,  the  integral  equation  suggests 
an  approximation  scheme  based  on  iteration  (Born  approximation).  In  general,  this 
should  succeed  if  the  scattered  wave  is  small,  as,  for  example,  when ‘the  wavelength 

*  With  our  choice  for  the  time  dependence,  (exp  +  ikr-iwt) /r  represents 
outgoing  (incoming)  spherical  waves. 
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is  large  compares  with  the  dimensions  of  the  scattered. 

As  a  simple  example  we  calculate  the  radiation  field  of  the  scattered  wave 
from  a  rigid  sphere  in  this  approximation  with  ine  aSSuapoidn  it*  «  i  ,  a-  using 
the  radius  of  the  sphere.  After  one  iteration  the  scattered  wave  is  given  approx¬ 


imately  by 


■Q  A** 

ps  ~ 


II  in=  ‘Pi 


*I«J 

(r  )  - - ] 

0  3ro  |r-t  |  V* 

o 


where  p_^  is  the  incident  wave,  which  we  assume  to  be  plane 


There  are  two  approximations.  We  keep  only  the--  terms  (the  radiation  field) 
and  we  expand  in  powers  of  ka  keeping  only  the  first  non-zero  contribution. 


■s  4ir 


dO  A  e 
o 


4-  4  4  **  4  . 

ik*r  „  exp  ik(r-r»r  +  ...) 

O  o  O 

3r  r  +  ... 


Aa  a 
4tt 


ikr  tt  2  2 

■ -  dn  {1  +  ika  cosy  +  ...H-ika  cos<fr  -  (ka)  cos  4>+. 

r  JJ  o 


where 


cos(r,ro) 


cos(k,ro) 


Carets  indicate  unit  vectors.  To  perform  the  integral  we  substitute  y  ■  0— , 


where  0  is  the  angle  of  r  measured  from  the  direction  of  the  incident  wave,  k  . 
The  first  non- vanishing  integral  gives 


a  ikr 
~  A £  e___ 

^  3  r 


(ka;  {l-cos8} 


-  57  - 


A  typical  diftmilerj  prs>M$»  is  that  gt  a  wav*  j#*§l8g  tg?eugh  m  t&erti 8« 
in  an  infinite  screen.  Eq.  (9)  is  usually  used  for  the  field  behind  the  screen 
together  with  the  Kirchhoff  approximation  vhich  assumes  that  the  values  »f-?w 
and  3pu/3n  on  the  aperture  are  just  what  would  be  given  by  the  incident  wave 
if  the  screen  were  not  present.  Then  the  integral  equation  is  well-defined. 
Methods  of  solution  are  discussed  in  most  books  on  acoustics,  optics,  and 
electromagnetic  theory. 

Appendix  I  :  The  Delta  Function 

The  delta  function  is  sometimes  "defined"  by  the  relations 

<5(x)  -  0  x  0 

b 

6(x)  dx  ■  1  a,b  >0  ^ 

-a 

from  which  immediately  follows  the  so-called  sifting  property: 

b 

[  f (x) 6 (x-x  )dx  -  f(x  )  a  <  x  <  b 

J  0  o 

a 

Intuitively,  it  is  a  function  which  is  zero  everywhere  except  at  it  =  0  where  it 
is  infinite  in  such  a  way  that  its  integral  is  unity.  It  is  clear  that  within  the 
context  of  ordinary  real  variable  theory  there  is  no  such  function,  for  a  function 
is  defined  only  if  it  is  assigned  a  definite  value  for  every  point  within  its 
range  (infinity  is  not  a  definite  value).  Moreover,  the  area  under  a  point  is 
always  zero. 

However,  tne  behavior  described  by  eqns.  (l)  can  be  approximated  with 
arbitrary  accuracy  by  ordinary  functions.  For  example,  consider  the  sequence 
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The  area  under  f  is  always  unity  and  the  maximum  at  x  =  0  "becomes  arbitrarily 
sharp  as  X  is  made  large. 

Therefore,  we  let  5 (x)  represent  the  limiting  behavior  of  such  a  sequence  and 
write,  symbolically 

5(x)  »  lim  f,(x) 

X-*»  A 

Identities  involving  derivatives  and  integrals  of  &  (x),  e«g.,  the  sifting 
property,  can  be  proved  provided  that  the  interchange  of  these  operations  with 
the  limit  can  be  Justified.  Justification  can  only  be  found  within  a  broader 
mathematical  context  than  the  ordinary  theory  of  real  variables.  Several  such 
contexts  are  available,  of  which  the  best  known  is  probably  the  theory  of  distribu¬ 
tions  as  presented  by  Lighthill  (1958),  whose  book  can  be  read  with  a  knowledge 
of  calculus  only.  For  an  intuitive  approach  we  also  recommend  the  book  of  van  detr 
Pol  and  Bremmer  (1964,  chap.  5).  Our  purpose  in  this  appendix  is  not  to  supply 
■rigor,  but  only  to  discuss  the  formal  manipulation  of  delta  functions. 

The  derivative  of  f  is  a  function  having 

A 

two  sharp  and  arbitrarily  close  extrema  as  x 
is  made  large.  It  too  haB  a  sifting  property. 

It  picks  the  values  of  a  function  at  two  points  which  are  arbitrarily  close  and, 
therefore,  produces  the  derivative.  By  a  formal  integration  by  parts  we  have 
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[  £(x}ba'*(x-xM)  dx  »  -f’{x) 

j  ~r-^'  o  o 

'  and,  generally,  _ 

j  f(x)  '«  (-l)n  £(n)  (xq) 

-  The  following  are- useful  identities; 

6(-x)  “  <5(x) 

5  (ax).  -  6(x) 

x  6(x)  «  0 


x  5’ (x)  *  ~d(x) 


-6(x?ra?)  -  {6  (x- a)  +  5(xta)} 


du 

dx 


d(x) 


U(x) 


x  >  0 
x  <  0 


By  the  equality  in  these  relations  we  mean  that  the  right  and  left  hand  sides  give 
the  same  result  when  substituted  into  an  integral.  •  Thus,  to  "prove"  the  second 
relation  we  form  the  integral  of  6  (ax)  with  an  arbitrary  continuous  function, 
f(x)  say,  and  ahov,  by  a  simple  change  of  the  variable  of  integration,  that  the 
result  is  |a|  *  f(0)  ,  which  is  what  is  given  by  the  right  hand  side  of  the 

identitiy. 

This  interpretation  of  the  equality  is  sufficient  for  our  purposes,  since  our 
delta  functions  always  ultimately  disappear  by  being  substituted  into  integrals  and, 
thus,  do  not  appear  in  the  final  results  of  our  calculations.  Their  use  can  be 
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regarded  as  an  intermediate  shortcut.  For  example,  eq.  <1.4)  can  he  derived  without 
the  use  of  the  delta  function  as  follows,  let  G  be.  a  harmonic  function  in  a 
region  v  except,  for  a  singularity  of  the  type  |r-ro|’^  at  one  point. 

Green’s  theorem  caxi  be  applied  to  the  region  v  -  vQ  which  excludes  from  v  a  small 
sphere  vQ  of  radius  B0  centered  at  the  singular  point  r  . 

[  dTo{G7Wo)  -  |  di0  <G  |f-  -  »  ||-} 

o  o 

v-v  0+0 

O  0 

In  the  limit  R  +  0  we  have 
o 


lia  |  R^  d O  '  {G  ~ 
R-H)J  °  0  9o 

o  ao  - 


-4tt  #(r) 


where  only  the  second  term  has  contributed,  and  the  result  is  established  by 
o  2 

substituting  Vq<2>  «  -4ir  p(rQ) ,  V*G  ■  0  in  the  volume  integral  over  v-vq.  ' 


Fourier  and  Laplace  Transforms 


The  Fourier  integral  theorem 


00  00 

f(x)  -  27  j  dp  |  dt  f(t) 

—00  —00 


may  he  written  in  two  steps: 


F(p) 

f(x) 


(1) 


(2) 


_1 

2it 


—00 
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(3) 


These  fora  a  transform  pair.  The  function  ?  (p)  is  called  the  Fourier  transform  of 
f  (x),  and  the  operator  in  eq.  (3)  through  -which  the  original  function  is  recovered 
Is  called  the 'inverse  Fourier  transform.  The  following  notation  is  also  used: 

F(p)  #  F{f>  =  f? 

*■*% -  }■  -  f(x) 

In  forming  the  transform,  pair  the  position  of  the  factor  of  2 it  and  of  the  plus 
and  minus  signs  in  the  exponentials  is  arbitrary.  Since  all  possible  combinations 
are  in  use  in  the  literature,  definitions  should  be  carefully  checked  when  using 
tables  of  Fourier  transforms. 

The  classical  theorem  for  the  validity  of  eq.  (l)  is  the'  Plsheherel 

theorem  which  states,  essentially,  that  a  sufficient  condition  is  the  existence  of 

& 

dxjf (x) ) ■  ,,  or,,  more  precisely,  that  f  (x)  is  Lebesgue  square  integrable 

4  —00 


(f  e  L  (-,-)) 


.  However,  this  theorem  is  too  restrictive  if  we  wish  to 


admit  generalized  functions,  such  as  the  delta  function  and  its  derivatives. 
Therefore,  we  shall  use,  instead,  a  sufficient  condition  given  by  Lighthill  (1958, 


p.2l)>  viz.,  the  existence  of  [  dx 

(14xZr 

This  allows  functions  bounded  by  polynomials  as  )x|  « 


for  some  integer  N. 


For  example, 


F”1 


F  I  •'1Plt  «»*  -  F 


F{1}  ■  2tt  6(p) 


Integrating  by  parts  we  have 


F  A  {6*  (p)}  ■  1x/2tt 
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or 


F{x)  *  -2iri  $’(p) 


arid,  in  general, 

F{xn}.  -  2ir(-i)a  n-  integer  >_  0 

Perhaps  the  most  important  property  of  the  Fourier  transform  for 
is  the  way  it  operates  on  derivatives.  Operating  on  eq.  (3)  with  the 
we  have 

f  *  *  F  *  {-ijjfp} 

or 

F{f '}  - 

and  hy  repetition 

F{f"}  -  -p2f  etc. 

P 

Thus,  a  differential  equation  is  transformed  to  an  algebraic  one. 

These  and  some  additional  properties  are  summarized  in  the  following  table. 


our  purposes 
derivative 
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IS- 


■[ 

l 


i 

| 

i 

i 

i 

* 

j 

i 

i 

i 

j 


t 

t 


! 

) 


»  «* 


f(x)  -  ^  |  dp  .-1I“  fp 

t  »  j  dx  eipx  f(x) 

— <30 

—as 

xn  n  ■  integer  ^  0 

2ir(-i)Vn)  (p) 

ix  f<x) 

^  fp 

dxp 

e  0  f (x) 

fp+P0 

-ipx 

f (x+XQ) 

e  fp 

f  (axtb) 

1  -ipb/a  . 

nr 8  v. 

f(x)  g(x) 

n  £p  *  bp 

f(x)  *  g(x> 

fpS 

The  star  indicates  the  convolution  of  two  function,  which  is  defined  by 


f(x)*g(x) 


09 

f(5)  g(x-5)  d5  »  g  *  f 

i 

—  00 


Our  sufficient  condition  does  not  admit  functions  which  are  exponentially 
large  as  x  -*  «  .We  can  allow  such  functions,  and  in  the  process 

define  the  Laplace  transform,  if  we  are  willing  to  throw  away  half  of  the  function 
and  assume  it  to  vanish  for  x  <  0  .  Thus,  we  assume 


5 
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f(x) 


■  0 


x  <  0 


I  f  (x)  I  <  eax  x>x>0  a>0 

0 

£he  Fourier  transform  of  f  (x)  may  not  exist,  we  define  a  new  function 

4>(x)  =  f(x)e  where  b  >  a 


«5fete6i-  In  exponentially  small  as-  -x  -i  If  we-. first  re~define  the 

transform  pair  by*  interchanging* /the* .-signs-'.in*  the*  exponentials  of 
<2,  3)  we  may  write 


hat  a  *•  b  +  ip 
txanaiota  of  $  . 


F{$} 


[  e  ipt  ij>(t)dt 


|  .-(w-tp)t  £(t)dt 
0 


.  As  a  function  of  p  F{$}  is  the  Fourier 

As  a  function* of  .s  we  define  irto'be  the- Laplace  transform  of  f 


f(t)dt 


We  obtain  the  inversion  theorem  through  the  inversion  theorem  for  the 
Fourier  transform  according  to  which 


$U)  "  f(x)e 
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-bx 
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2tt 


dP  .lp:t  ♦ 
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b+i" 


dp  e(b+lp>\ 
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2tt1 


ds  es 


b-i“ 


there  was  originally  no  upper  limit  for  b  so  that  the  contour  for  the  inversion 
integral, C  say,  which  is  a  straight  line  parallel  to  the  imaginary  axis  of  the 
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s-plane,  is  arbitrarily  far  to  the  right  of  the  origin,  but  by  Cauchy's  theorem 
is  equivalent  to  any  straight  line  to  the  right  of  all  singularities.  For  x  >  0 
the  integral  can  be  evaluated  in  terms  of  contributions  from  the  poles,  branch 
cuts,  etc.  to  the  left  of  C  .  For  x  <  0  the  inversion  integral  automatically 
gives  zero  since  the  integrand  vanishes  as  C  is  moved  to  the  right. 

Integration  by  parts  gives  the  effect  of  the  Laplace  transform  on 
derivatives,  which  we  include  in  the  following  table. 

f(x)  |  =  fs 


For  the  Laplace  transform  convolution  is  defined  by 

x 

f  *  g  -  j  d£  f«>  g<x-5)  -  g  *  f 
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SOLUTIONS 


Chapter  II 


Suppose  {G(r,r^)  -  G(r,ro)}ds  f  0.  Since  G  is  determined  to  within 

an  additive  constant,  let  F(r^)  «  j  ds  G(r,r^)  and  replace  GCr,^)  by 

G(r,rp  -  F(r^)/s  where  s  is  the  total  area. 
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Chap tar  IV 
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